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ABSTRACT 


Radiation  dose  distribution  measurements  have  been  carried  out  by  the  Naval  Surface 
Warfare  Center  (NSWC)  using  200  and  391  MeV  electron  beams  completely  penetrating 
layers  of  Aluminum,  lead  and  aluminum,  polymethyl  methacrylate  (PMMA),  and  PMMA 
sandwiching  an  air  gap.  For  the  case  of  391  MeV  electrons,  the  results  have  been  compared 
previously  to  the  corresponding  distributions  predicted  by  the  Monte  Carlo  simulation  codes 
EGS3  and  ACCEPT  fRef  1].  Those  measurements/predictions  for  391  MeV  are  here 
compared  to  predictions  by  the  CYLTRAN  electron/photon  transport  code,  and  the  200  MeV 
measurements  done  by  NSWC  are  compared  to  CYLTRAN  calculations.  The  CYLTRAN 
code  predictions  agree  well  with  measurements  at  391  MeV.  Comparison  of  CYLTRAN 
calculations  with  the  NSWC  results  for  200  MeV  indicates  possible  saturation  of  the 
detectors  used  to  take  the  measurements.  The  distribution  of  energy  dose  within  the  target 
has  a  large  dependence  on  the  location  of  the  air  gap.  The  variation  in  dose  distribution  is 
caused  by  the  change  in  target  geometry  resulting  from  insertion  of  the  air  gap. 
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I .  INTRODUCTION 

The  objective  of  this  work  is  to  test  the  ability  of  a 
computer  model,  using  the  Integrated  Tiger  Series  (ITS)  code, 
to  accurately  predict  the  behavior  of  radiation  cascades 
caused  when  a  high  energy  electron  beam  traverses  a  solid 
target.  The  approach  used  has  been  to  study  experimental 
measurements,  construct  a  computer  model  suitable  for  use  as 
input  to  the  ITS  simulation  code  system,  run  the  code,  and 
compare  the  results  to  the  measurements.  In  order  to 
understand  the  approach  taken  in  the  construction  of  the 
computer  model  forming  the  basis  of  this  work,  some  details  of 
the  experiments  to  which  the  model  will  be  compared  must  be 
examined . 

A.  BACKGROUND 

Researchers  from  the  Naval  Surface  Warfare  Center  (NSWC) 
have  conducted  experiments  at  the  Massachusetts  Institute  of 
Technology,  Bates  Linear  Accelerator  Center  to  measure  the 
three-dimensional  spatial  distribution  of  energy  deposited  by 
radiation  showers  in  a  multilayer,  multi-element  target  stack 
when  the  target  is  bombarded  by  a  high  energy  electron  beam 
[Ref.  1,  2,  and  3].  The  experiment  by  NSWC  was  conducted  using 
200  MeV  and  391  MeV  electrons  incident  on  four  separate  target 
stack  configurations.  The  results  of  the  391  MeV  experiment 
were  compared  with  energy  deposition  distributions  calculated 
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using  two  different  Monte  Carlo  electron  transport  model 
computer  codes,  EGS3  (Electron-Gamma  Shower,  Version  3)  and 
ACCEPT.  The  EGS  (for  Electron-Gamma-Shower)  system  of 
computer  codes  is  a  general  purpose  package  developed  by  the 
Stanford  Linear  Accelerator  Center  (SLAC)  for  the  Monte  Carlo 
simulation  of  coupled  electron  and  photon  transport  in 
arbitrary  geometry.  ACCEPT  is  one  of  a  series  of  codes  based 
on  ETRAN  [Ref.  5]  called  the  Integrated  Tiger  Series  (ITS) , 
developed  by  Sandia  National  Laboratory  [Ref.  4]. 

The  target  materials  examined  in  the  NSWC  wor]:  were 
Aluminum,  Lead  followed  by  Aluminum,  Polymethyl  Methacrylate 
(PMMA,  also  known  as  Lucite  or  Plexiglass) ,  and  PMMA 
sandwiching  an  air  gap.  Target  plates  of  roughly  equal 
thickness  were  constructed  of  each  type  of  target  material. 
These  5  or  6  plates  were  stacked  together  to  make  target 
stacks.  Each  target  stack  was  configured  so  that  the  total 
target  thickness,  expressed  as  the  areal  mass  density,  was 
approximately  50  g/cm^.  A  schematic  of  a  target  stack  as 
configured  for  the  NSWC  experiments  is  shown  in  Figure  1. 

Interspaced  at  regular  intervals  in  the  target  stack  were 
planar  arrays  of  Lithium  Fluoride  Thermoluminescent  Dosimeters 
(TLDs) .  The  TLD  arrays  consisted  of  88  individual  dosimeters, 
0.1cm  in  diameter,  0.3cm  in  length,  arranged  in  concentric 
annular  patterns  with  their  centers  along  the  beam  axis,  and 
were  designed  so  as  to  measure  the  amount  of  energy  deposited 
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(dose)  as  a  function  of  radial  distance  from  the  beam  at 
several  positions  along  the  length  of  the  target  stack. 

The  target  stacks  were  irradiated  by  electron  beams  of  200 
and  391  MeV.  The  results  obtained  in  the  work  by  the  NSWC  form 
the  basis  of  this  work  and  are  reported  in  references  1  and  2 
for  391  MeV,  and  in  Reference  3  for  200  MeV. 

The  work  by  NSWC  included  a  comparison  of  the  391  MeV 
measurements  to  predictions  by  two  Monte  Carlo  electron/photon 
shower  simulation  cedes,  EGS3  and  ITS/ACCEPT.  This  comparison 
led  the  res.-iarchers  at  NSWC  to  three  specific  conclusions 
(these  conclusions  pertained  only  to  the  391  MeV  data) .  The 
NSWC  conclusions  are  summarized  below.  [Ref.  1-3] 

•  CYLTRAN  and  EGS3  both  predicted  values  for  the  dose 
distribution  that  agreed  well  with  measurements  except  in 
the  PMMA  target  configurations  with  an  air  gap. 

•  Neither  ACCEPT  nor  EGS3  was  able  to  accurately  predict  the 
dose  distribution  in  the  PMMA  configurations  with  an  air 
gap.  The  codes  predicted  doses  as  much  as  100%  larger  than 
the  measurements. 

•  The  dose  distribution  in  the  PMMA  target  stacks  was 
largely  dependent  on  the  location  of  the  air  gap  within 
the  stack. 

In  comparison  of  the  results  from  the  present  work  to 
those  from  the  NSWC  work,  the  following  conclusions  have  been 
reached. 

•  CYLTRAN  predictions  agreed  well  with  the  391  MeV  data 
taken  by  NSWC  with  no  exceptions.  Even  though  the  dose 
distribution  in  the  target  configuration  with  an  air  gap 
varied  with  the  location  of  the  air  gap,  CYLTRAN  was  able 
to  predict  the  dose. 

•  The  variation  in  dose  distribution  caused  by  changing  the 
location  of  the  air  gap  seems  to  be  a  geometry  effect. 
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This  conclusion  was  drawn  from  the  results  obtained  by 
comparing  the  dose  in  a  target  with  an  air  gap  to  that 
obtained  when  the  air  gap  was  replaced  by  a  void,  and 
finding  that  they  were  identical. 

•  A  comparison  of  the  NSWC  200  MeV  data  to  CYLTRAN 
calculations  was  done.  This  was  a  new  comparison  since  the 
NSWC  work  had  not  compared  the  data  to  a  transport  code  at 
this  energy.  CYLTRAN  produced  results  that  agree  very  well 
with  the  data  at  200  MeV  except  for  on-axis  measurements 
in  two  target  configurations  (PMMA  and  Lead  followed  by 
Aluminum) .  Analysis  of  these  results  indicates  that  the 
dosimeters  used  to  take  the  data  may  have  been  saturated. 


B.  SOME  CONSIDERATIONS  IN  ELECTRON  TRANSPORT  MODELING 
1.  Electron  Step  Size 

The  number  of  collisions  undergone  by  an  electron  as 
it  passes  through  a  target  are  too  numerous,  except  possibly 
at  very  low  electron  energies,  to  allow  treatment  of 
individual  collisions  even  by  the  most  powerful  computers. 
Instead,  the  electron  path  through  the  material  is  divided 
into  discreet  segments  in  which  a  large  number  of  collisions 
occur.  The  smaller  pathlengths,  or  steps,  allow  the 
calculation  of  energy  deposition  by  the  sampling  of  net 
angular  deflection  and  net  energy  loss  using  the  appropriate 
multiple  scattering  approximation  (the  Goudsmit-Saunderson 
approximation  for  the  ITS  CYLTRAN  code  [Ref.  5])  and  a 
"condensed  history  model"  approach  [Ref.  6]. 

Seltzer  explains  that  the  step  size  is  chosen  by  the 
code  user  by  weighing  a  set  of  conflicting  requirements.  The 
steps  should  be  short  enough  that  (1)  most  of  the  steps  in  the 
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electron  history  are  completely  inside  the  target  so  that  the 
use  of  multiple-scattering  theories  for  unbounded  media  are 
appropriate,  (2)  the  average  energy  loss  is  small  so  that  the 
use  of  one-velocity  multiple-scattering  distributions  is 
justified,  and  (3)  the  net  angular  deflection  is,  on  the 
average,  small  so  that  the  spatial  displacement  for  the  step 
can  be  approximated  closely  by  a  straight  line.  The  steps 
should  be  long  enough  that  (1)  there  are  a  sufficient  number 
of  collisions  per  step  to  justify  the  use  of  multiple¬ 
scattering  theories,  and  (2)  the  number  of  steps  per  history 
is  not  too  large.  [Ref.  7] 

The  portion  of  the  ITS  code  package  that  generates 
cross-section  data  for  the  transport  code  to  use  in  its 
calcula<-ions  is  called  XGEN.  The  step  size  is  determined 
during  the  execution  of  XGEN,  in  the  cross-section  generation 
process,  so  that  it  satisfies  the  relation 

The  value  k=8  is  chosen  by  default  in  XGEN,  the  ITS 
cross-section  generating  program,  so  that  in  passing  through 
a  given  step,  an  electron  will  lose  8.3%  of  its  energy  [Ref. 
8]  . 

Translation  of  the  percentage  decrease  in  electron 
energy  to  an  equivalent  physical  step  length  is  done  by 
assuming  the  electron  slows  down  continuously  over  the  step. 
The  step  size  is  given  as 
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(2) 


where  is  the  electron  range  at  the  initial  and  final 
energies  in  the  step  as  given  by  the  continuous-slowing-down 
approximation . 

2.  Electron  Radiation  Length 

It  is  convenient,  in  considering  the  passage  of 
electromagnetic  radiation  through  material,  to  measure  the 
thickness  of  the  material  in  terms  of  the  radiation  length.  In 
most  processes,  some  or  all  of  the  dependence  on  the  medium  is 
contained  in  the  radiation  length.  The  radiation  length  is 
defined  as  that  distance  Lp  over  which  a  high  energy  electron 
is  slowed  to  a  fraction  equal  to  1/e  of  its  original  energy. 
The  values  of  radiation  length  for  electrons  for  materials  of 
interest  in  this  work  are  given  in  Table  I.  The  table  gives 
values  of  radiation  length  in  units  of  cm  and  g/cm^,  the 
latter  being  the  product  of  the  former  and  the  density  of  the 
material . 

Table  I.  RADIATION  LENGTHS  FOR  SELECTED  MATERIALS  [Ref.  9] 


Material  Rad.  Length  (cm)  Rad.  Length 


Aluminum  8.9  24.01 
Lead  0.56  6.37 
PMMA  34.4  40.55 
Air  30420  36.66 


3.  Thin  Film  Scattering 

The  experimental  work  on  which  this  work  is  based 
involved  a  target  stack  configuration  with  an  air  gap  between 
two  sections  of  solid  material  (Lucite) .  The  possibility  was 
explored  in  the  present  work  that  the  ITS  code  might  better 
predict  doses  in  air  if  the  air  gap  were  treated  as  a  thin 
film,  due  to  the  low  density  of  air,  and  the  long  radiation 
length  of  an  electron  in  air.  Most  electron/photon  transport 
codes  in  current  use,  including  the  ITS  series  of  codes,  are 
subject  to  the  constraint  that  they  must  be  used  in  situations 
where  the  material  depth  exceeds  the  step  size  calculated  for 
the  electron  as  it  traverses  the  material.  The  Goudsmit- 
Saunderson  theory  breaks  down  in  the  limit  of  small  material 
thicknesses  (thin  films) .  Work  in  this  area  has  produced  the 
Jordan-Mack  correction  to  the  CYLTRAN  code  to  correct  the 
multiple  scattering  distribution  in  the  code  to  give  better 
agreement  with  measurements  for  thin  film  targets.  [Ref. 3,  pp. 
14-16]  The  effect  of  the  Jordan-Mack  correction,  which  can  be 
applied  to  the  CYLTRAN  code  by  the  user,  was  investigated  in 
this  case  to  determine  if  the  same  sort  of  thin  film 
phenomenon  might  be  involved  in  electron  transport  through  an 
air  gap. 

C.  OVERVIEW  OF  ITS 

The  Integrated  Tiger  Series  of  computer  codes  was 
developed  over  the  last  20  years  at  the  U.S.  Department  of 
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Energy's  Sandia  National  Laboratory  for  use  in  the  simulation 
of  electron/photon  transport  through  various  materials.  ITS  is 
based  primarily  on  ETRAN  (for  Electron  TRANSport) ,  which  was 
developed  by  the  National  Bureau  of  Standards  [Ref.  7,  p. 
249].  The  codes  forming  the  basis  of  ITS  are  TIGER,  CYLTRAN, 
and  ACCEPT.  The  three  codes  differ  primarily  in  that  (1)  they 
treat  problems  in  different  dimensionalities  and  (2)  they  are 
designed  to  handle  different  problem  geometries.  The  TIGER 
code  is  confined  to  the  analysis  of  problems  in  one  dimension. 
The  three-dimensional  codes,  CYLTRAN  and  ACCEPT,  have  been 
applied  in  this  work  and  will  be  explained  in  further  detail 
below. 

1 .  CYLTRAN 

CYLTRAN,  a  fully  three-dimensional  transport  code 
employing  cylindrical  geometry,  is  a  natural  choice  for 
simulating  problems  involving  electron  beam  sources.  A 
complete  description  of  the  ITS-CYLTRAN  code  system  is 
provided  in  Appendix  A. 

CYLTRAN  enables  the  code  user,  relatively  easily,  to 
break  a  target  with  cylindrically  symmetric  geometry  into  an 
arbitrary  number  of  zones  and  will  determine  the  energy 
deposited  in  each  zone  by  a  beam  of  electrons  (as  in  this 
case)  or  photons.  The  decision  of  the  fineness  of  division  of 
the  target  into  zones  must  be  balanced  against  the  need  for 
timely  results,  as  the  cost  for  more  (or  finer)  zones  is 
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longer  run  time  in  executing  the  CYLTRAN  code.  The  CYLTRAN 
geometry  allows  the  electron  beam  to  be  incident  on  the 
cylinder  at  any  arbitrary  point  and  from  any  direction.  The 
CYLTRAN  code  has  no  ability  to  present  results  in  a  graphical 
form.  The  user  must  provide  the  resources  to  take  the  output 
from  CYLTRAN  (which  by  default  is  manageable,  but  can  be  quite 
lengthy)  and  manipulate  it  either  manually  or  electronically 
into  a  form  which  can  be  analyzed  graphically  or  in  some  other 
manner. 

In  this  work,  the  CYLTRAN  code  was  given  an  input  file 
which  gave  geometry  data  defining  the  boundaries  of  the 
desired  number  of  target  zones.  The  ensuing  code  calculations 
produced  a  formatted  output  file  containing  the  energy 
deposition  distributions  and  other  code  results.  This  output 
file  was  read  by  an  interactive  program  written  in  the 
Interactive  Data  Language  (IDL) ,  which  manipulated  the 
appropriate  geometry  data  and  energy  deposition  distributions 
into  plottable  formats  for  analysis. 

2 .  ACCEPT 

ACCEPT  is  a  fully  three-dimensional  electron/photon 
transport  code  which  is  applicable  to  problems  of  any  geometry 
using  a  combinatorial  geometry  scheme.  While  the  combinatorial 
geometry  system  makes  ACCEPT  more  complex  to  use,  its  basic 
computational  algorithms  are  identical  to  those  in  CYLTRAN  and 
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the  two  should  give  identical  results,  as  long  as  the  problem 
geometry  allows  the  use  of  either  code. 
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II.  SIMULATION  PROCEDURE 


A.  GENERAL  PROCEDURAL  APPROACH 

As  mentioned  previously,  references  1  through  3  reported 
the  results  of  measurements  taken  by  NSWC,  and  also  reported 
the  results  of  the  comparison  of  measurements  at  391  MeV  beam 
energy  with  dosages  predicted  by  two  different  Monte  Carlo 
electron-photon  transport  simulation  codes,  EGS3  and  ACCEPT. 
The  measurements  taken  by  NSWC,  reported  in  reference  3,  were 
not  compared  to  any  computer  code  calculations. 

The  NSWC  results  reported  in  references  1  through  3  were 
modeled  in  this  work  using  the  CYLTRAN  member  of  the  ITS  code 
package.  The  CYLTRAN  code  was  run  on  a  DEC  VAXStation  3200 
minicomputer  using  VAX  Fortran.  The  results  of  the  CYLTRAN 
runs  were  read  and  plotted  by  a  program  written  for  this  work 
using  the  Interactive  Data  Language  (IDL)  ,  a  commercially 
available  graphics  package  marketed  by  Research  Systems, 
Incorporated.  A  user's  guide  and  the  source  code  for  the  IDL 
program  written  for  the  manipulation  and  plotting  of  the 
calculations  done  in  this  work  are  attached  as  Appendix  B. 

Although  the  target  stacks  used  in  the  NSWC  experiment 
were  actually  rectangular,  or  box-shaped,  the  homogeneous 
nature  of  the  target  stack  materials  and  the  use  of  a  nonnally 
incident,  monoenergetic  electron  beam  suggested  that  the 
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targets  be  modeled  as  a  right  circular  cylinder  with  the 
electron  beam  (which  is  taken  to  be  a  "point"  beam)  traveling 
along  the  cylinder's  longitudinal  axis  and  incident  on  the 
face  of  the  cylinder.  This  geometry  allowed  the  natural  use 
of  the  CYLTRAN  member  of  the  ITS  code  with  its  simpler  input 
stream.  Figure  2  is  a  schematic  of  the  target  model  in  the 
case  of  the  NSWC  target  consisting  of  PMJ-IA  sandwiching  an  air 
gap  located  downstream  of  the  fi'-st  10  g/cm^  of  target 
material . 

The  first  runs  of  the  CYLTRAN  code  were  intended  to 
determine  whether  CYLTRAN  in  fact  produced  the  same  results, 
within  statistical  expectations,  as  those  obtained  at  NSWC 
using  ACCEPT.  The  onief  concern  in  these  runs  was  the 
prediction  of  the  dose  in  the  TLDs  interleaved  between  each 
target  'slab'  along  the  length  of  the  target,  as  opposed  to 
later  simulation  runs,  which  would  divide  the  target  stack 
into  many  more  zones,  attempting  to  resolve  in  greater  detail 
what  was  happening  to  the  electron  energies  as  the  particles 
progress  through  the  target.  For  reasons  of  distinguishing 
between  the  degree  of  fineness  of  the  zones,  these  firsc 
target  geometries  were  called  "low  resolution  targets"  and  the 
corresponding  CYLTRAN  runs  were  called  "low  resolution"  runs. 

After  the  initial  low  resolution  CYLTRAN  runs  had  been 
completed,  a  series  of  runs  was  done  to  determine  how  the  dose 
in  the  air  gap  and  the  layers  of  PMMA  surrounding  it  responded 
to  changes  in  various  model  characteristics.  The 
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Figure  2.  Low  Resolution  Target  With  a  39cin  Air  Gap  Downstream 
of  One  PMMA  Targe  c  Section.  This  target  is  called  a  low 
resolution  target  because  tie  target  material  itself  is  not 
resolved  into  zones  in  which  the  dose  is  evaluated. 


character!  .cs  which  were  varied  for  comparison  were  the 
number  of  CYLTRAN  input  zones,  the  number  of  primary  electron 
histories  to  be  followed  by  CYLTRAN,  the  scattering  angle 
approximation  used  in  the  code,  and  the  location  of  the  air 
gap. 

CYLTRAN  allows  the  user  to  determine  the  number  and 
location  of  spatial  zones  to  be  monitored  within  the  target 
cylinder.  The  user  simply  places  into  the  input  data  file 
geometry  data  specifying  the  location  of  zone  boundaries  (in 
the  terminology  used  by  the  ITS  system,  each  zone  is  called  an 
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"input  zone") .  One  can  treat  the  entire  target  cylinder  as  a 
single  input  zone  or  divide  the  target  cylinder  into  hundreds 
of  finely  divided  input  zones  and  track  the  dose  deposited  in 
each  zones,  if  so  desired.  In  this  series  of  runs,  the  CYLTRAN 
model  cylinder  and  its  target  sections  were  divided  into 
small,  highly  resolved  "input  zones".  Each  section  of  target 
material  was  divided  longitudinally  into  10  slices  of  equal 
thickness,  and  the  39cm  air  gap  was  divided  into  13  slices, 
each  having  a  thickness  of  3cm.  Each  slice  in  all  materials 
was  itself  divided  into  12  radial  zones  to  identically  match 
the  radial  locations  of  the  LiF  TLD's  used  in  the  NSWC 
experiments.  The  resultant  large  number  of  input  zones  (828  in 
the  case  of  PMMA  sandwiching  a  39cm  air  gap,  660  for  all  other 
configurations  considered  here)  allowed  CYLTRAN  to  track  the 
energy  deposition  at  a  much  finer  scale,  and  a  "high 
resolution"  view  was  obtained  of  the  energy  dose  in  the  target 
as  a  function  of  radius  and  depth.  Figures  3  and  4  compare  the 
zoning  of  the  low  resolution  and  high  resolution  cases. 

The  high  resolution  zoning  approach  to  the  CYLTRAN  model 
was  applied  to  all  targets  of  concern  in  the  NSWC  experiment 
in  order  to  better  assess  the  performance  of  CYLTRAN  in 
predicting  the  measured  dose. 

Some  CYLTRAN  runs  were  done  to  test  the  plausibility  of 
the  idea  that  the  air  gap  could  be  treated  as  a  thin  film.  The 
Jordan-Mack  thin  film  scattering  correction  to  the  ITS  codes 
was  applied  to  CYLTRAN.  Runs  were  conducted  on  one  of  the 
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et  Section  and  Dosimeter  Array 


targets  of  concern.  The  results  of  this  run  were  compared  to 
the  results  for  the  same  target  obtained  from  CYLTRAN  without 
the  correction  applied  to  assess  the  correction's  effect,  if 
any,  on  the  dose  distribution  in  the  air  gap. 

It  was  of  some  interest  in  this  work  to  evaluate  the 
statistical  uncertainty  involved  in  the  CYLTRAN  calculations, 
and  to  get  some  idea  of  how  this  uncertainty  changed  with  the 
choice  of  input  parameters.  Of  particular  interest  was  the 
question  of  how  the  CYLTRAN-calculated  dose  in  the  target 
varies  with  a  change  in  the  number  of  primary  electrons.  To 
this  end,  it  was  necessary  to  consider  how  many  electrons 
would  be  sufficient  to  statistically  represent  the  real 
experiment.  In  the  NSWC  experiments,  the  targets  were  exposed 
for  10  to  20  seconds  to  a  O.l^A  current.  This  translates  to  a 
total  of  from  1.6x10’^  to  3.2xl0’^  electrons.  Diagnostic 
CYLTRAN  runs  on  the  Digital  VaxStation  3200  used  in  this  work 
indicated  that,  with  a  baseline  run  time  of  about  2  hours  for 
10,000  primary  electrons,  increasing  the  number  of  primary 
electrons  10  fold  to  100,000  would  require  a  run  time  of  about 
24  hours,  or  a  factor  of  10  increase  in  computer  run  time. 
Changing  the  number  of  primary  electrons  from  10,000  to 
100,000  causes  a  decrease  in  uncertainty  from  20-50%  in  the 
former  case  to  10-20%  in  the  latter  (for  all  materials  in  this 
work  except  air) .  Due  to  relatively  low  value  of  uncertainty 
for  100,000  electrons,  and  to  the  fact  that  any  further  gains 
in  the  level  of  certainty  would  cost  dearly  in  computer  run 
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time,  a  CYLTRAN  run  with  100,000  electrons  was  decided  to  be 
sufficiently  representative  of  the  real  beam.  The  CYLTRAN  runs 
for  the  high  resolution  target  were  repeated  using  100,000 
electron  histories  and  the  impact  of  the  gain  in  statistical 
certainty  was  assessed. 

B.  NEW  ANALYSIS  FOR  200  MEV  ELECTRONS 

Reference  3  listed  the  measurements  obtained  in 
experiments  using  200  MeV  electron  beams.  Data  was  reported  by 
NSWC  on  the  three  target  stacks  consisting  of  only  Aluminum, 
Lead  followed  by  Aluminum,  and  PMMA.  Computer  simulations  had 
not  been  done  previously  against  this  data.  As  part  of  this 
work,  CYLTRAN  runs  were  performed  to  see  how  well  CYLTRAN 
would  predict  the  dose  distribution  at  this  energy. 
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III.  RESULTS  AND  DISCUSSIONS  OF  PREVIOUSLY  ANALYZED  DATA  AT 
391  MEV 

A.  RESULTS  WITH  CYLTRAN  CORRECTED  FOR  THIN  FILM  SCATTERING 

The  Jordan-Mack  correction  for  use  in  correcting  the  ITS 
code  to  accurately  predict  energy  deposition  in  thin  films, 
when  applied  to  CYLTRAN  in  this  work,  has  no  effect  on  the 
dose  distribution  in  any  of  the  target  configurations 
considered.  This  is  an  indication  that,  if  the  codes  are 
inherently  unable  to  treat  an  air  gap  at  high  energies,  that 
inability  is  not  related  to  any  inability  to  handle  thin  film 
scattering  phenomena. 

B.  GENERAL  COMPARISON  OF  CYLTRAN  TO  ACCEPT  AND  NSWC  DATA 

Reference  1  reported  dose  distributions  in  the  target 
stack  as  predicted  by  ACCEPT.  ACCEPT  calculations  were 
reported  for  dose  in  the  LiF  dosimeters  at  6  longitudinal 
points  along  the  stack  in  the  case  of  the  2  targets  with  PMMA 
sandwiching  an  air  gap,  and  for  5  points  along  the  stack  in 
the  cases  of  the  other  targets  (PMMA,  aluminum,  and  lead 
followed  by  aluminum) .  Figures  5  through  9  show  the  results  of 
predictions  calculated  using  CYLTRAN  (it  is  important  to  note 
that  figures  5  through  9  show  points  calculated  at  a  few 
discreet  points  along  the  length  of  the  cylinder  corresponding 
to  the  longitudinal  depths  at  which  the  NSWC  data  were  taken. 
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Lines  drawn  between  the  points  are  shown  only  to  help 
distinguish  between  sets  of  calculations/data) .  These 
calculations  were  done  using  the  "low  resolution”  model 
detailed  in  Chapter  II  so  that  a  fair  comparison  could  be  done 
between  the  two  sets  of  calculations.  Generally,  agreement 
between  the  NSWC  data  and  CYLTRAN  calculations  is  good  to  very 
good.  Discrepancies  between  the  code  and  experiment  seems  to 
be  greatest  at  a  radius  of  5  cm.  In  view  of  a  statistical 
uncertainty  at  this  radius  of  20  -  40  percent,  however,  the 
discrepancy  may  not  be  significant.  ITS  calculations,  in  all 
target  materials,  agree  very  well  for  radii  in  the  range  from 
the  axis  to  about  2  cm. 

CYLTRAN  scores  a  dose  by  whether  it  was  deposited  by  a 
primary,  knock-on,  or  gamma-secondary  electron.  Analysis  of 
the  CYLTRAN  calculations  shows  that  dose  on  axis  is  due  almost 
totally  to  primary  electrons.  At  a  radius  of  20  cm,  the  dose 
is  due  mostly  to  gamma-secondary  electrons.  The  doses  in 
between  are  due  to  various  combinations  of  primary  and  gamma¬ 
secondary  electrons.  Knock-on  electrons  have  only  a  minor 
contribution  in  all  cases. 

ACCEPT  and  CYLTRAN  give  results  that  vary  by  little  in 
most  cases.  However,  in  some  cases,  especially  as  the  radius 
approaches  20cm,  the  two  codes  seem  to  give  widely  varying 
answers.  It  has  been  noted  previously  that  ACCEPT  and  CYLTRAN 
have  the  same  transport  physics  as  their  basis,  so  their  is  no 
reason  to  believe  the  two  codes  should  produce  widely  varying 
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Figure  5.  Low  Resolution  391  MeV  CYLTRAN  Dose  Distribution  in 
PMMA  with  39cin  Air  Gap  Downstream  of  1st  PMMA  Target  Section, 
10,000  primary  electrons.  Note  close  agreement  of  CYLTRAN  and 
data.  Compare  with  Figure  10. 
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Figure  6.  Low  Resolution  391  MeV  CYLTRAN  Dose  Distribution  in 
PMMA  with  39cin  Air  Gap  Downstream  of  4th  PMMA  Target  Section, 
10,000  primary  electrons. 
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Figure  9.  Low  Resolution  391  MeV  CYLTRAN  Dose  Distribution  in 
Target  of  Lead  followed  by  Aluminum,  10,000  primary  electrons. 
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results.  However,  when  consideration  is  given  to  the  fact 
that,  at  20cin  radius  the  uncertainty  in  the  calculations  can 
be  as  high  as  50%,  the  variation  between  the  results  is  not 
that  surprising.  That  is,  ACCEPT  and  CYLTRAN  seem  to  be  within 
the  statistical  uncertainty  of  each  other. 

C.  EFFECT  OF  AIR  GAP  LOCATION  ON  TARGET  DOSE  DISTRIBUTION 

If  an  air  gap  is  inserted  in  the  target  stack,  there  is  a 
marked  variation  in  the  dose  distribution.  Th^  dose  in  the 
portion  of  the  target  material  on  or  close  to  the  beam  axis 
doesn't  change  significantly.  However,  as  one  looks  at  the 
results  at  distances  away  from  the  axis,  the  dose  in  any 
particular  section  of  target  material  is  significantly 
affected  by  whether  or  not  it  is  preceded  in  the  longitudinal 
direction  by  an  air  gap.  As  the  air  gap  is  moved  to  a 
different  depth  in  the  target  stack,  the  dose  in  the  material 
preceding  the  air  gap  is  reduced.  However,  the  section  of 
target  material  immediately  following  the  air  gap  receives  a 
higher  dose  than  that  which  would  have  been  received  if  the 
air  gap  were  not  there.  Both  of  these  effects  are  increasingly 
pronounced  as  radius  from  the  beam  axis  increases.  A 
comparison  of  Figures  5,  6,  and  7  illustrates  this  point. 

CYLTRAN  calculations  agree  very  well  with  data  in 
reference  1.  That  is,  the  changes  in  dose  due  to  insertion  of 
an  air  gap  seem  to  be  faithfully  predicted  within  statistical 
uncertainty,  which  will  be  discussed  later.  Incidentally,  in 
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those  cases  where  there  are  noticeable  discrepancies  between 
CYLTRAN  and  data,  the  CYLTRAN  dose  distribution  mows  a 
tendency  to  overpredict.  The  only  case  in  which  the  agreement 
between  CYLTRAN  and  the  data  may  be  questionable  is  shown  in 
Figure  5  at  5cm  radius  and,  most  remarkably,  at  20cm  radius. 
In  this  regard,  it  is  well  to  note  the  large  statistical 
uncertainty  inherent  in  the  CYLTRAN  Monte  Carlo  algorithm, 
which  runs  from  25-50  percent  at  20cm  radius  for  10,000 
primary  electrons  in  the  target  of  Figure  5,  as  will  be 
discussed  later.  In  light  of  this  uncertainty,  the  apparer  : 
axsagreement  in  Figure  5  between  CYLTRAN  and  data  may  not  be 
significant . 

The  low  resolution  CYLTRAN  ca?culations  snown  in  figures 
5  through  9  show  only  the  doses  in  the  dosimeter  arrays 
located  between  the  target  sections  and  before  and  after  the 
air  gap.  A  high  resolution  CYLTRAN  run,  in  which  target 
sections  are  divided  in  much  finer  zones,  sheds  a  great  deal 
of  light  on  how  the  dose  distribution  behaves  in  the  target 
material  itself.  Figure  10  illustrates  the  high  resolution 
dose  distribution  in  the  case  of  PMMA  with  the  air  gap  located 
downstream  of  the  first  target  section.  Figure  11  shows  the 
dose  distribution  for  PMMA  with  the  air  gap  downstream  of  the 
fourth  target  section.  Figures  12,  13,  and  1^'  show  dose 
distributions  for  targets  of  solid  PMMA,  Aluminum,  and  Lead 
followed  by  Aluminum,  respectively. 
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A  comparison  of  Figure  10  with  Figure  5  shows  that  the 
overall  disagreement  between  the  NSWC  data  and  these  CYLTRAN 
calculations  is  reduced  by  dividing  the  problem  target 
geometry  into  smaller  zones.  The  high  resolution  calculations 
verify  the  continuous  nature  of  the  dose  distribution.  Except 
for  statistical  fluctuations  of  the  dose  about  a  mean  value, 
the  dose  is  in  fact  smoothly  distributed  in  the  cylinder,  even 
across  the  air  gap.  The  most  well-defined  results  are  given  on 
axis.  This  is  what  would  be  expected,  as  there  the  beam  is 
most  intense  and  the  probability  of  an  electron  colliding  with 
an  "air  molecule"  would  be  the  highest.  In  other  words,  the 
statistics  are  better  on  axis  and  get  increasingly  worse  with 
increasing  radius.  This  would  explain  the  absence  of  plottable 
points  at  5cm  and  20cm  radius  in  the  air  gap  in  Figure  10.  At 
radii  greater  than  about  5cm  in  the  air  gap  with  10000  primary 
electrons  the  probability  of  a  collision  is  extremely  low. 
Thus,  the  dose  in  the  air  gap  will  be  calculated  by  CYLTRAN  to 
be  zero  in  most  cases.  This  is  born  out  in  Figure  10. 

D.  NUMBER  OF  PRIMARY  ELECTRON  HISTORIES 

Since  the  Monte  Carlo  calculations  which  are  the  basis  of 
the  CYLTRAN  code  is  statistical  in  nature,  one  might  expect 
that  increasing  the  number  of  primary  electrons  would  cause  a 
decrease  in  the  statistical  uncertainty  of  the  CYLTRAN 
calculations.  Results  obtained  in  this  work  would  indicate 
that  this  is  so.  A  plot  of  dose  per  incident  electron  for  the 
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target  of  Figure  10  using  100,000  incident  electrons  is  shown 
as  Figure  15.  The  better  definition  of  the  dose  distribution 
at  large  radii  in  the  air  gap  is  a  certain  indication  that 
statistics  have  improved  significantly.  Comparison  of  figures 
10  and  15  reveals  that,  in  the  target  material,  the  value  of 
the  dose  calculated  by  CYLTRAN  doesn't  depend  on  the  number  of 
primary  electron  histories  followed.  The  statistical  benefit 
in  raising  the  number  of  histories  is  shown  in  the  comparison 
of  Figures  16  and  17.  A  factor  of  10  increase  in  the  n,’mber  of 
primary  electrons  produces  a  decrease  in  uncertainty  of  at 
least  50%  at  all  radii.  The  most  significant  benefit  derived 
by  using  more  electrons  is  gained  in  the  air  gap  and  at  large 
radii,  where  the  collision  probability  is  smallest. 

E.  DOSE  BREAKDOWN  BY  ELECTRON  TYPE 

The  dose  deposited  in  a  target  zone  by  an  electron  can  be 
classified  by  the  type  of  electron  that  deposited  the  dose. 
The  electron  type  is  identified  as  either  primary,  gamma¬ 
secondary,  or  knock-on.  A  primary  electron  is  an  electron  that 
was  in  the  original  electron  beam  incident  on  the  target.  A 
knock-on  electron  is  one  that  has  been  knocked  free  of  an  atom 
through  some  collisional  process.  A  gamma-secondary  electron 
is  an  electron  that  has  been  emitted  from  an  atom  through  any 
process  involving  a  photon.  An  important  part  of  this 
classification  is  that  once  a  photon  has  knocked  an  electron 
from  an  atom,  that  electron  is  a  gamma-secondary  electron  as 
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Figure  10.  High  Resolution  391  MeV  CYLTRAN  Dose  Distribution 
in  PMMA  with  39cin  Air  Gap  Downstream  of  1st  PMMA  Target 
Section,  10,000  primary  electrons.  Compare  to  figure  5.  Note 
reduced  disagreement  with  data. 
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Figure  11.  High  Resolution  391  MeV  CYLTRAN  Dose  Distribution 
for  PMMA  with  39cin  Air  Gap  Downstream  of  4th  PMMA  Target 
Section,  10,000  Primary  Electrons.  CYLTRAN  agrees  very  well 
with  data  here. 
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Figure  12.  High  Resolution  391  MeV  CYLTRAN  Dose  Distribution 
for  PMMA  Target,  10,000  Primary  Electrons.  CYLTRAN  agrees  very 
well  with  data. 
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Figure  14.  High  Resolution  391  MeV  CYLTRAN  Dose  Distribution 
for  Target  of  Lead  Followed  by  Aluminum,  10,000  Primary 
Electrons.  CYLTRAN  agrees  very  well  with  data. 
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Figure  15.  High  Res.  391  MeV  CYLTRAN  Dose  Distribution  in  PMMA 
with  Air  Gap  Downstream  of  1st  PMMA  Target  Section,  100,000 
primary  electrons.  Compare  to  previous  figure  to  note  drastic 
decrease  in  statistical  uncertainty. 
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Figure  16.  Uncertainty  in  Total  Dose  by  Radius  for  10,000 
Primary  Electrons  in  Target  of  Figure  9.  Compare  to  following 
figure  where  100,000  electrons  are  used.  Note  the  improvement 
in  statistics  with  a  large  number  of  electrons. 
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Figure  17.  Uncertainty  in  Total  Dose  by  Radius  for  100,000 
Primary  Electrons  in  Target  of  Figure  10.  Compare  with  the 
previous  figure  for  10,000  electrons.  Note  the  improvement  in 
statistics  with  a  larger  number  of  primary  electrons. 
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well  as  all  electrons  on  down  the  chain  that  the  secondary 
electron  might  start.  Figure  18  illustrates  the  classification 
of  these  dose  deposition  mechanisms. 

The  CYLTRAN  output  provides  detailed  analysis  of  whether 
the  dose  in  any  particular  target  zone  comes  from  a  primary, 
gamma-secondary,  or  knock-on  electron.  Figures  19  through  23 
show  the  dose  per  incident  primary  electron  at  various  radii 
on  the  target  of  Figure  10  with  100,000  primary  electrons. 
Figures  19  and  2  0  show  that,  on  axis,  the  dose  is  almost 
totally  due  to  loss  of  energy  by  primary  electrons.  The 
negative  dose  contributed  by  the  knock-on  electrons  simply 
says  that  a  small  portion  of  the  energy  deposited  by  the 
primaries  is  carried  out  of  the  zone  by  knock-ons  created  in 
the  zone.  Note  that  by  the  time  the  primary  electrons  traverse 
the  first  PMMA  target  section  and  enter  the  air  gap,  the  on- 
axis  dose  is  reduced  to  25%  of  its  original  value.  Once  the 
electrons  are  one  third  of  the  way  through  the  air  gap,  the 
on-axis  dose  due  to  primary  electrons  has  been  reduced  to 
about  5%  of  its  original  value.  Figure  20  shows  that,  although 
the  dose  on-axis  through  the  end  of  the  air  gap  is  mostly  from 
primaries,  immediately  downstream  of  the  air  gap  the  dose 
becomes  dominated  by  gamma-secondary  electrons.  The  small 
energy  contribution  from  knock-ons  does  not  play  a  significant 
role  on-axis  except  in  the  first  target  section.  We  see  from 
the  off-axis  CYLTRAN  calculations  in  Figures  21  and  22  that, 
while  there  is  very  little  dose  away  from  the  axis  in  the 
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Figure  18.  Classification  of  Dose  Deposition  by  Type  of 
Electron  Depositing  the  Dose. 

first  target  section,  as  the  beam  progresses  down  the  cylinder 
the  contribution  from  the  primaries  is  Increasingly  overtaken 
by  that  from  gamma-secondaries  farther  away  from  the  beam 
axis.  Finally,  as  shown  in  Figure  23,  the  dose  at  a  radius  of 
20cm  is  almost  totally  due  to  gamma-secondary  electrons.  In 
fact,  there  is  very  little  energy  deposited  in  the  cylinder  at 
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Figure  19.  Dose  Breakdown  by  Electron  Type  on  Beam  Axis  for 
Target  of  Figure  10  with  100,000  Electrons,  391  MeV.  Most  of 
the  dose  is  due  to  primary  electrons  in  the  first  15cm  depth. 
The  next  figure  expands  the  region  beyond  2Ccm  depth. 
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Figure  20.  Dose  Breakdown  on  Beam  Axis  (Expanded  View  Beyond 
20cm  Depth)  for  100,000  Electrons  at  391  MeV.  Photon  processes 
cause  rapicJ  increase  in  gamma-secondaries  at  air  gap  boundary, 
so  beyond  48cm  dose  is  mostly  from  low  energy  g-secondaries. 


41 


c  I39|e ;  b)(  /  Aen )  esoQ 


Figure  21.  Dose  Breakdown  by  Electron  Type  at  2cin  Radius  for 
Target  of  Figure  10  with  100,000  Electrons,  391  MeV.  Dose  at 
this  radius  is  mostly  from  primaries  and  is  highest  at  mid¬ 
depth  after  electrons  have  scattered  to  that  radius. 
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Figure  22.  Dose  Breakdown  by  Electron  Type  at  Scm  Radius  for 
Target  of  Figure  10  with  100,000  Electrons,  391  MeV. 
Statistics  in  the  air  gap  begin  to  decrease  because  few 
electrons  reach  that  radius  inside  the  gap. 
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Figure  23.  Dose  Breakdown  by  Electron  Type  at  20cin  Radius  for 
Target  of  Figure  10  with  100,000  Electrons,  391  MeV.  Dose  from 
knock-ons  overshadows  others  in  the  air  gap.  Uncertainty  in 
the  air  gap  is  very  high  due  to  the  few  collisions  occurring. 


44 


20cm  until  the  air  gap  is  passed,  at  which  point  the  dose  due 
to  secondary  and  knock-on  electrons  rises  dramatically,  the 
knock-on  dose  falling  slowly  back  to  a  very  low  level  and  the 
secondary  dose  falling  and  becoming  a  constant  at  a  level 
slightly  higher  than  the  knock-on  dose.  This  dose  breakdown 
leads  to  the  conclusion  that  the  rise  one  sees  in  the  dose 
after  the  air  gap  as  indicated  in  Figure  10  is  formed  by  an 
increase  contribution  from  gamma-secondary  electrons  somehow 
built  up  in  the  beam's  passage  through  the  air  gap. 

Calculations  identical  to  those  above  for  the  case  of 
10,000  electrons  have  been  done,  the  results  of  which  are 
attached  in  Appendix  C  as  Figures  Cl  through  C4 . 

F.  EFFECT  OF  REPLACING  THE  AIR  GAP  BY  A  VOID 

Some  insight  into  the  nature  of  the  calculated  dose 
distribution  in  the  air  gap  was  gained  by  substituting  a  void 
for  the  39cm  air  gap.  The  result  of  doing  this  is  shown  as 
Figure  24.  The  dose  in  the  void  is,  of  course,  zero.  But, 
perhaps  more  significantly,  there  is  no  appreciable  change  in 
the  dose  distribution  in  the  target  material  in  this  case  from 
that  reflected  in  Figure  10,  the  identical  target  with  an  air 
gap.  These  results  would  tend  to  suggest  that  the  increased 
dose  in  the  target  material  downstream  of  the  air  gap  is  not 
an  artifact  of  the  air  itself,  but  rather  one  of  the  target 
geometry . 
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Figure  24.  CYLTRAN  Dose  Distribution  for  Target  of  Figure  10 
with  a  39cm  Void  Replacing  the  Air  Gap.  This  figure  and  figure 
10  are  very  similar,  showing  that  dose  distribution  dependence 
on  air  gap  location  is  a  geometry  effect. 
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6.  AIR  GAP  STEP  SIZE  EFFECTS 


Work  done  by  Bielajew  and  Rogers  on  calculations  of  the 
dose  distributions  in  air  due  to  high  energy  electrons  has 
shown  that  if  the  default  Monte  Carlo  step  size  is  used,  one 
could  expect  the  dose  to  be  overpredicted  [Ref.  10].  In  this 
work,  it  has  not  been  possible  to  investigate  the  effect  of 
changing  the  step  size,  as  CYLTRAN  does  not  have  the 
capability  of  accomplishing  this  for  only  one  material  in  a 
multi-material  target.  However,  at  an  energy  of  391  MeV,  the 
step  size  required  to  obtain  the  optimum  results,  according  to 
Biejalew  and  Rogers,  would  be  about  one  target  diameter,  or 
20cm.  This  step  size  would  require  the  value  of  k  in  equation 
1  to  be  equal  to  about  19,400. 

H.  RANDOM  NUMBER  EFFECTS 

The  dose  in  the  material  has  a  large  dependence  on  the 
value  of  the  initial  random  number  seed  used  to  start  the  ITS 
calculation.  Obvious  though  this  observation  may  be  for  a 
Monte  Carlo-based  algorithm,  it  is  well  to  keep  in  mind  that 
any  result  achieved  using  the  code  is  subject  to  a  calculable 
uncertainty.  Figures  C5  through  C7  in  Appendix  C  are  given  as 
examples  of  the  effect  of  changing  the  random  number  seeds. 
When  compared  to  their  counterparts,  they  reflect  noticeable 
variations.  However,  the  variations  are  verifiably  within  that 
statistical  uncertainty  which  is  given  by  CYLTRAN  along  with 
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the  dose  values.  The  uncertainty  in  a  dose  in  solid  materials 
in  this  work  typically  runs  about  20-50%. 
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IV.  RESULTS  FOR  PREVIOUSLY  UNRMRLYZED  200  MEV  DATA 


CYLTRAN  calculations  were  done  to  check  the  validity  of 


predictions 

at  200 

MeV. 

Calculations  were  compared 

to 

measurements 

made  at 

Bates  and  reported 

in  Reference 

1. 

Figures  25, 

26,  and 

27 

illustrate  the 

results  of 

the 

comparison  for  Aluminum,  PMMA,  and  Lead  followed  by  Aluminum 
targets,  respectively. 

The  results  of  CYLTRAN  runs  compare  very  well  with  the 
measurements  in  all  three  targets,  as  was  expected  from  the 
previous  results  (discussed  in  Chapter  III)  for  the  391  MeV 
data.  The  sole  notable  deviation  of  CYLTRAN  from  the  data  is 
in  the  latter  depths  of  the  Pb/Al  and  PMMA  targets.  In  both  of 
these  cases,  the  CYLTRAN  calculations  on  the  beam  axis  begin 
to  provide  results  lower  than  the  data  in  the  third  target 
section.  The  deviation  on  axis  increases  continuously  through 
the  end  of  the  problem  cylinder.  Comparison  of  figures  25,  26, 
and  27  to  figures  12,  13,  and  14  in  chapter  III  (which  are  the 
corresponding  problems  at  391  MeV) ,  shows  that  general  trends 
in  the  CYLTRAN  calculations  for  the  391  MeV  case  have  been 
repeated  in  the  200  MeV  case,  as  one  might  expect.  Also, 
examination  of  the  NSWC  data  as  shown  in  figures  25  and  27 
(particularly  for  Pb/Al  and  PMMA  targets)  shows  that  the  on- 
axis  dose  data  begins  to  plateau  after  about  two  target 
sections  depth.  This  seems  suspect  since,  in  all  cases  looked 
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at  in  this  work,  after  a  dose  peak  in  the  first  few 
centimeters  of  target  material,  dose  generally  decreases 
continuously.  Along  with  those  facts,  it  is  noteworthy  that 
the  CYLTRAN  calculations  have  produced  very  good  agreement  in 
every  other  case  examined  in  this  work.  There  is  no  reason  to 
expect  such  a  significant  difference  in  the  results  for  the 
two  energies,  since  the  target  configurations  involved  are 
identical.  All  this  would  tend  to  suggest  that  the  TLD's  used 
in  the  NSWC  measurements  for  the  target  configurations  of  PMMA 
and  Lead  followed  by  Aluminum,  on  axis  at  least,  had  become 
saturated  and  did  not  provide  accurate  data. 
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Figure  25.  High  Res.  CYLTRAN  Dose  Dist.  for  200  MeV  Electrons 
in  PMMA  Target.  CYLTRAN  agrees  well  with  data  except  on  beam 
axis  and  large  depth-  Linearity  of  NSWC  data  past  25cm  depth 
suggests  saturation  of  TLD*s  on  axis.  Compare  with  figure  12. 
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Figure  26.  High  Resolution  CYLTRAN  Dose  Distribution  for  200 
MeV  Electrons  in  Aluminum  Target,  10,000  Primary  Electrons. 
CYLTRAN  agrees  well  with  data  for  this  case.  Compare  to  figure 
13  at  391  MeV. 


Figure  27.  High  Res.  CYLTRAN  Dose  Distribution  for  10,000  200 
MeV  Electrons  in  Target  of  Lead  followed  by  Aluminum.  CYLTRAN 
agrees  with  data  except  on  axis  at  large  depth.  Again,  TLD 
saturation  is  suggested.  Compare  to  figure  14. 
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V.  CONCIUSIONS 


In  general  the  results  obtained  at  Bates,  as  detailed  in 
Reference  1,  have  been  verified.  A  comparison  of  the  Bates 
data  to  the  calculations  done  here  using  CYLTRAN  yields  the 
following  conclusions. 


•  The  dose  in  any  particular  section  of  target  material  is 
significantly  affected  by  whether  or  not  it  is  preceded  in 
the  longitudinal  direction  by  an  air  gap.  However,  this 
dependence  is  not  a  product  of  the  air  itself.  Rather,  it 
is  an  effect  of  the  change  in  geometry  produced  when  the 
air  gap  is  inserted.  We  come  to  this  conclusion  by 
replacing  the  air  gap  by  a  void  and  getting  identical 
results.  A  small  angular  deflection  gives  a  small  lateral 
displacement  at  short  range,  but  gives  a  much  larger 
lateral  displacement  at  longer  ranges.  See  Chapter  III, 
Section  C. 

•  Uncertainty  in  dose  calculations  in  the  air  gap  is  very 
large  due  to  the  small  number  of  collisions  undergone  by 
an  electron  as  it  traverses  the  air  gap.  The  uncertainty 
is  as  much  as  99%  in  some  cases.  This  is  not  affected  by 
changing  the  number  of  input  zones  (dividing  the  air  more 
finely) .  However,  it  is  decreased  significantly  by 
increasing  the  number  of  primary  electrons  used  in  the 
simulation.  See  Chapter  III,  Section  D. 

•  In  the  PMMA/Air/PMMA  target  stack,  the  dose  on  axis  is  due 
almost  totally  to  primary  electrons.  At  a  radius  of  20  cm, 
the  dose  is  due  mostly  to  gamma-secondary  electrons.  The 
doser  in  between  are  due  to  various  combinations  of 
primary  and  gamma-secondary  electrons.  Knock-on  electrons 
have  only  a  minor  contribution  in  all  cases.  See  Chapter 
III,  Section  E. 

•  In  constructing  a  model,  the  numerical  results  are  not 
affected  by  the  degree  of  fineness  with  which  the  target 
is  divided  into  input  zones.  However,  the  added  sense  of 
continuity  in  the  data  provided  when  the  problem  geometry 
is  more  finely  divided  are  tremendously  beneficial.  See 
Chapter  III,  Section  C. 
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•  Correcting  the  ITS  code  to  treat  the  air  gap  as  a  thin 
film  seems  to  have  no  affect  on  the  dose  distribution.  See 
Chapter  III,  Section  A. 

•  The  final  value  of  dose  in  the  material  depends  greatly  on 
the  value  of  the  initial  random  number  seed  used  to  start 
the  ITS  calculation.  That  is,  any  one  ITS  calculation  is 
not  absolute  in  its  correctness.  The  value  given  for  the 
energy  deposited  in  a  particular  zone  in  this  work  could 
vary  from  20-50%.  However,  the  results  obtained  using 
different  random  seeds  are  statistically  consistent,  given 
that  the  statistical  uncertaii,ty  in  the  calculations  is 
also  20-50%.  See  Chapter  III,  Section  H. 

•  CYLTRAN  agrees  well  with  data  for  all  of  the  targets 
considered  in  this  work.  This  is  true  for  electron 
energies  at  both  200  MeV  and  391  MeV.  The  391  MeV  data  has 
been  analyzed  before  by  NSWC  using  ACCEPT.  Our  analysis 
shows  the  failure  of  the  previous  ACCEPT  comparison  to  the 
data  is  probably  due  to  lack  of  resolution  and  a  small 
number  of  electron  histories. 

•  The  analysis  of  the  200  MeV  data  has  not  been  done 
previously.  This  work  shows  that  CYLTRAN  is  able  to 
accurately  calculate  the  measured  dose.  Beyond  a  radiation 
length,  the  dose  on  axis  begins  to  deviate  from  the  data. 
The  trends  in  the  data  compared  to  the  calculations 
suggest  that  the  TLD's  used  in  taking  the  200  MeV 
measurements  were  saturated. 
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APPENDIX  A.  USING  ITS  ON  THE  DEC  VAXStation  3200  COMPUTER 


In  this  work,  ITS  was  installed  and  executed  on  a 
dedicated  DEC  VAXStation  3200  computer.  The  work  has  not  been 
completed  without  a  significant  amount  of  time  spent  in  self- 
education  in  the  DEC  terminology  and  the  use  of  the  VAX  VMS 
operating  system  used  by  the  computer  itself,  and  in  learning 
the  subtleties  of  properly  using  the  ITS  computer  code.  In  the 
hope  of  saving  future  ITS  code  users  some  time  at  the  front  of 
the  "learning  curve",  this  appendix  will  attempt  to  explain 
rudimentary  details  of  how  to  (1)  turn  on  (boot)  the  computer, 
log  on  to  a  user  account  and  to  begin  to  work,  and  (2)  how  to 
execute  a  CYLTRAN  run  using  the  ITS  system  as  currently 
installed  on  the  VAXStation  3200. 

A.  BOOTING  AND  USING  THE  VAXSTATION  3200 
1.  Turning  the  Computer  On 

Turn  on  the  Constant  Voltage  Conditioner  with  the  red 
power  switch.  It  will  make  an  annoying  constant  humming  noise, 
which  is  nothing  to  be  alarmed  about. 

Turn  on  the  computer  by  pressing  the  switch  on  the 
right  of  the  little  VAXStation  3200  logo  on  the  front  panel  of 
the  computer  cabinet.  The  switch  will  glow  orange  after  being 
pressed. 
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The  system  will  go  through  diagnostics  and  booting 
procedures  for  a  few  minutes.  The  entire  booting  procedure 
should  be  complete  in  less  than  five  minutes. 

After  the  computer  has  completed  booting,  the  monitor 
screen  should  show  the  DECWindows  login  screen  with 
distinctive  huge  blocked  letters  spelling  out  DIGITAL  and  a 
small  window  below  containing  spaces  to  type  in  the  username 
and  password.  At  this  time  users  are  able  to  log  on  to  the 
system,  if  a  user  account  has  been  set  up  by  accessing  the 
computer's  System  Manager  account.  If  a  user  account  has  not 
been  previously  established,  it  is  recommended  that  it  be  done 
immediately  by  accessing  the  System  Manager  account  as 
explained  in  Section  2  below  and  by  following  instructions  in 
the  computer's  system  management  manual  (a  set  of  several 
volumes,  usually  kept  on  shelves  next  to  the  computer 
station) .  It  is  a  dangerous  habit  (speaking  from  personal 
experience)  to  use  the  System  Manager  account  for  routine  user 
tasks . 

2 .  Logging  on  to  the  Computer 

After  the  computer  has  completed  the  booting  process, 
you  are  ready  to  log  on  to  the  computer  and  begin  work.  First, 
press  any  button.  If  it  wasn't  awake  already,  the  monitor  will 
"wake  up"  and  display  the  DECWindows  login  screen  (described 
in  Section  1) .  Select  USER  by  clicking  the  left  mouse  button 
on  the  space  for  the  username  and  type  your  username.  Then 
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select  PASSWORD  and  type  your  password.  Then  press  the  ENTER 
key.  The  computer  will  begin  to  log  in  to  your  account  and 
will  put  the  DECWindows  user  screen  on  the  monitor.  This 
process  takes  a  couple  of  minutes.  After  you  log  in,  you  can 
enter  commands  to  the  computer  by  using  the  pull-down  menus 
(not  a  very  fast  way  to  do  it)  or  by  clicking  on  the  DECterm 
window  (which  activates  the  window  so  you  can  enter  a  command) 
and  type  commands  directly  (much  better) . 

If  you  wish  to  log  in  to  SYSTEM  (the  System  Manager 
account) ,  type  the  username  SYSTEM  and  the  password  ESMERALDA. 
You  must  log  on  to  SYSTEM  to  establish  any  necessary  user 
accounts.  There  will  be  times,  since  this  is  a  machine 
essentially  dedicated  to  one  user,  that  the  user  will  have  to 
log  on  as  SYSTEM  to  accomplish  some  task  that  can  only  be  done 
by  the  System  Manager,  such  as  device  allocation  and  memory 
management.  After  you  log  on  to  SYSTEM  and  the  work  screen 
appears  on  the  monitor,  click  on  the  DECterm  window  to  enter 
commands  at  the  "$"  prompt. 

To  log  out,  select  SESSION  on  the  Session  Manager 
Window  Menu  Bar  and  click  on  QUIT.  You  will  be  asked  to 
confirm  that  you  want  to  quit  and,  if  so,  you  will  be  logged 
out . 

3.  Turning  the  Computer  Off 

Log  in  to  SYSTEM  as  described  above.  Click  on  the 
DECterminal  window  and  type  SHUTDOWN. 
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When  the  system  goes  completely  through  the  shutdown 
sequence,  which  takes  a  minute  or  two,  the  message  SYSTEM 
SHUTDOWN  COMPLETE — USE  CONSOLE  TO  HALT  SYSTEM  appears  on  the 
monitor.  Shut  off  the  computer  with  the  orange  switch  and  then 
turn  off  the  Constant  Voltage  Conditioner  with  the  red  power 
switch. 

4.  Miscellaneous  Comments  on  Operation  of  the  Computer 

Users  accomplish  tasks  on  the  computer  through  the 
entering  of  user  commands.  Commands  are  entered  to  the 
VAXStation  in  a  DEC  language  called  DCL  (Digital  Command 
Language) .  Rather  than  try  to  explain  a  lot  of  commands  here, 
it's  best  that  you  review  the  VAXStation  User  Manual  volume 
titled  "DCL  Dictionary"  to  get  a  feel  for  the  commands  and 
what  they  can  do.  A  few  commands  that  apply  directly  to  the 
use  of  ITS  and  IDL  (Interactive  Data  Language)  will  be 
mentioned  in  this  appendix  and  in  Appendix  B. 

The  VAXStation  is  currently  connected  to  an  ethernet 
accessing  several  computers  in  the  Naval  Postgraduate  School 
Physics  Department  and  other  campus  facilities,  including  the 
Physics  Simulation  Lab  and  the  W.R.  Church  Computer  Center. 
Access  to  this  network  is  through  a  smaller  network  (the 
DECnet)  of  DEC  VAXes  (there  are  currently  4  VAXes  on  the 
DECnet)  belonging  to  the  Physics  Department.  To  access  another 
node  on  the  network,  the  DCL  command  is 
SET  HOST  nodename 
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You  must  know  the  nodename  and  it  must  be  defined  in 
the  local  database.  If  you  have  the  right  nodename  the  other 
computer  will  prompt  you  for  a  username  and  password. 

There  are  two  user-interface  systems  available  on  the 
VAXStation:  DECWindows  and  VWS  (VAX  Window  System) .  Currently 
the  default  on  this  machine  is  set  to  DECWindows  due  to  the 
user  friendliness  of  the  interface.  The  VWS  is  a  more 
barebones  window  system  that  runs  faster  but  doesn't  make 
things  especially  easy  for  the  inexperienced  user.  The  window 
system  can  be  changed  by  making  changes  in  the  SYSTEM  account 
file  SYSTARTUP_V5.COM  and  going  through  a  reboot  procedure. 
Read  the  manual  on  the  topic  or  seek  help  from  a  more 
experienced  user. 

Remember  to  do  complete  system  backups  on  a  regular 
basis.  To  backup  the  entire  system; 

•  Set  the  Break  Enable/Disable  switch  to  Enable  (in  back  of 
the  computer  above  the  LED  number  display) .  Point  the 
switch  toward  the  circle  with  a  dot  inside  it. 

•  Shutdown  the  system. 

•  When  you  get  the  message  SYSTEM  SHUTDOWN  COMPLETE — USE 
CONSOLE  TO  HALT  SYSTEM  appears  on  the  monitor  press  the 
halt  button  on  front  of  the  computer  twice.  You  should  get 
the  >>>  prompt. 

•  Type  the  command:  b/eOOOOOOO  (that's  seven  zeros). 

•  Write  protect  the  hard  disk  by  pressing  the  button  on 
front  of  the  computer. 

•  Perform  the  backup  with  the  command  (all  on  one  line) : 

"BACKUP/IMAGE/VERIFY/BUFFER_C0UNT=5  (space) 

DUAO ; MUAO : name . BCK/REWIND/LABEL=date 
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Load  a  blank  TK70  cartridge  in  the  tape  drive  when 
prompted.  When  the  backup  is  finished,  unload  the  tape, 
write  enable  the  hard  disk  and  restart  the  computer. 


B.  ITS  OVERVIEW 

The  Integrated  Tiger  Series  of  Electron/Photon  transport 
codes  (ITS)  is  the  most  widely  used  particle  transport  code  in 
the  world.  The  code  package  was  developed  to  incorporate  eight 
individual  codes  which  were  developed  over  the  period  from 
1968  to  1981.  All  the  codes  are  based  on  the  original  ETRAN 
model  developed  by  M.  Berger  and  S.  Seltzer.  The  ITS  code 
system  consists  of  four  primary  code  packages,  which  are 
listed  in  Table  Al. 

Table  Al.  FOUR  PRIMARY  ITS  CODE  PACKAGES. 

XDATA:  The  electron  and  photon  cross  section  file. 

XGEN:  The  cross  section  generation  program. 

ITS:  The  Monte  Carlo  program  file. 

UPEML:  A  machine  portable  update  emulator. 

The  heart  of  the  ITS  is  the  program  library  file  ITS, 
which  contains  the  eight  Monte  Carlo  programs  plus  system 
directives  for  the  CRAY,  IBM,  VAX,  and  CDC  operating  systems. 
The  update  emulator  program  UPEML  creates  the  various  Monte 
Carlo  codes  for  a  given  system  with  any  corrections  to  those 
codes  that  may  be  desired.  The  output  fortran  source  code  from 
UPEML  is  then  compiled,  linked,  and  stored  as  an  executable 
module.  Program  XGEN  generates  the  problem  specific  cross 


62 


section  data  tape  using  file  XDATA  for  referenced  inputs  and 
a  user  defined  input  file.  The  Monte  Carlo  codes  then  read  in 
the  cross  section  tape  and  process  the  user  defined  problem. 

One  of  the  eight  ITS  codes  is  CYLTRAN,  which  simulates  the 
transport  of  particle  trajectories  through  a  three-dimensional 
multimaterial  right  circular  cylinder.  For  this  project  only 
the  CYLTRAN  code  was  used.  As  an  ITS  user  the  following  steps 
were  required  to  execute  an  ITS  run. 

1.  Create  the  specific  ITS  code  CYLTRAN  with  the  required 
correction  schemes. 

2.  Generate  a  cross  section  tape  based  on  the  different  type 
of  materials  contained  in  the  cylindrical  geometry  of  a 
problem. 

3.  Create  an  input  file  which  lists  all  the  input  parameters 
required  to  calculate  desired  outputs. 

4.  Submit  the  input  file  and  the  generated  cross  section 
tape  to  the  ITS  Monte  Carlo  codes  to  execute  a  run. 

Table  A2  is  a  sample  input  file  to  generate  the  cross 
section  tape  for  the  materials  in  a  cylindrical  geometry.  Each 
material  line  represents  a  different  medium  in  the  cylinder. 
Percentages  of  each  material  in  a  compound  and  its  density 
must  be  specified.  Single  element  lines  such  as  Cu,  has  its 
density  stored  in  ITS  and  is  automatically  used  for  the 
simulation  when  needed. 

Once  a  cross  section  tape  is  generated,  an  input  file  with 
the  parameters  designed  for  a  particular  simulation  must  be 
created.  Section  D  below  contains  a  sample  input  file  to 
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Table  A2.  SAMPLE  INPUT  FILE  TO  CREATE  A  CROSS  SECTION  TAPE 


Energy  16.0 
Material  Au 
TITLE 

391  MeV  Cross  Section  for  Au  Foil 


execute  a  CYLTRAN  run  for  one  of  the  problems  in  this 
work. [Ref .  3,  pp.  39-40] 

C.  RUNNING  ITS  ON  THE  VAX8TATION  3200 

The  VAXStation  3200  is  currently  set  up  to  run  ITS  in 
either  a  batch  mode  or  interactively.  In  the  batch  mode,  the 
user  can  submit  the  code  run  to  the  computer  as  a  job  separate 
from  any  user  interactive  processes.  This  way,  the  user  can 
continue  to  use  the  computer  interactively  while  the  ITS  run 
is  being  processed  in  the  background  by  the  computer.  In  the 
interactive  mode,  the  user  cannot  use  the  computer  terminal 
while  the  run  is  in  progress.  The  batch  mode  is  by  far  the 
most  efficient  method  to  use  when  submitting  ITS  runs. 

Before  any  ITS  code  may  be  run  on  the  computer,  an 
executable  version  of  the  code  must  be  generated.  Assume  that 
you  have  generated,  using  the  UPEML  processor  provided  with 
the  ITS  package,  an  ITS  specific  fully  compilable  file  with 
the  name  [  .ITS]CYL0.FOR  which  you  wish  to  run  on  the  computer. 
Generate  the  executable  file  by  following  this  procedure: 


Enter  the  command  FOR  [ . ITS]CYL0. FOR.  The  computer's  VAX 
Fortran  compiler  will  compile  the  fortran  file  and  save  an 
object  file  in  your  default  directory  named  CYLO.OBJ. 


•  Enter  the  command  LINK  CYLO.OBJ.  The  computer  will  link 
the  program  properly  and  will  produce  an  executable  file 
named  CYLO.EXE  and  save  it  in  your  default  directory. 

•  The  executable  file  can  then  be  executed  under  one  of  the 
procedures  below. 

The  process  to  be  followed  in  submitting  an  ITS  run  is 
given  below.  The  example  given  is  typical  of  a  run  done  in 
this  work.  As  ITS  is  currently  configured  on  the  VAXStation 
3200,  the  process  must  be  followed  as  given.  In  the  example, 
it  is  assumed  that  a  specific  ITS  code  file,  a  cross  section 
file,  and  a  problem  input  file  have  been  generated. 

1.  Batch  Mode  Submission  of  an  ITS  Run 

Before  an  ITS  job  can  be  submitted  in  batch  mode,  two 
separate  data  files  must  be  constructed  using  the  VAX  text 
editor.  The  first  file  contains  1)  a  command  telling  the 
computer  where  to  look  for  input  required  in  processing  the 
batch  job,  and  2)  a  command  telling  the  computer  to  run  the 
executable  file  containing  the  ITS  code.  An  example  of  the 
first  file  is  the  file  BATCH.COM,  which  contains  the  following 
two  lines. 

DEFINE/USER_MODE  SYS$INPUT  ’ [ YAW] BATCH_IN . DAT ' 

RUN  [.ITSjCYLO_BATCH.EXE 


The  second  data  file  contains  the  name  of  the  XGEN- 
generated  cross  section  file  and  the  name  of  the  file 
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containing  the  problem  input  parameters  and  geometry  input.  An 
example  of  the  second  data  file  for  this  work  is  the  file 
BATCH_IN.DAT,  which  could  contain  the  following  two  lines. 

[ YAW . PMMA ] CYLl 0_3  9_4  0 . INP 
[ YAW . I TS  2 ] XT APE . OUT 

Once  the  two  required  files  have  been  created,  the 
batch  job  may  be  submitted.  Submit  the  job  by  entering  the  DCL 
command 

SUBMIT/NOTIFY/NOPRINTER  BATCH.COM 
where  the  filename  BATCH.COM  is  the  name  of  the  file 
containing  the  name  of  the  executable  ITS  file.  The  terms 
NOTIFY  and  NOPRINTER  are  keywords  telling  the  computer  to 
notify  you  when  the  batch  is  completed  and  to  send  the  job  log 
to  a  file  in  your  login  directory,  but  not  to  the  print  queue. 
See  the  DCL  Dictionary  for  other  keywords  to  use  with  the 
SUBMIT  command. 

After  submitting  the  batch  job,  you  may  return  to 
interactive  use  of  the  computer  while  the  computer  processes 
the  job  in  the  background.  The  nice  thing  about  the  batch 
process  is  that  you  may  queue  as  many  runs  as  you  want  and  let 
the  machine  run  overnight  or  over  a  weekend.  However,  before 
doing  this,  it's  a  smart  idea  to  set  the  batch  queue  job  limit 
to  1  by  logging  in  to  SYSTEM  and  entering  the  command 

SET  QUEUE/ J0B_LIMIT=1  SYS$BATCH 
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Be  sure  to  log  out  of  SYSTEM  and  back  in  to  your  user  account 
before  submitting  a  job.  If  the  job  limit  were  higher  than  1 
(say  3,  for  instance)  the  computer  would  try  to  run  all  3  jobs 
at  the  same  time,  and  the  run  time  for  all  3  jobs  would  be 
significantly  increased.  Also,  the  computer's  capabilities  for 
running  that  many  jobs  at  once  are  really  limited,  and  it 
could  bomb  out  altogether.  With  a  job  limit  of  1,  the  jobs  run 
one  at  a  time,  with  an  unlimited  number  of  jobs  held  in  the 
queue  pending  execution. 

2.  Running  ITS  Interactively  on  the  VAXStation  3200 

The  interactive  mode  is  the  most  user-friendly  way  to 
run  ITS  on  the  VAXStation  3200.  To  execute  the  ITS  run,  enter 
the  command  RUN  followed  by  the  filename  of  the  ITS  executable 
code  you  wish  to  run.  A  command  to  execute  a  run  in  this  work 
would  look  like  this: 

RUN  [ .ITS] CYL0.EXE 

where  the  file  [.  ITS] CYL0.EXE  is  an  executable  ITS  problem 
specific  code  generated  for  this  work. 

The  computer  will  then  begin  executing  the  code  and 
interactive  use  of  the  computer  will  no  longer  be  available. 
After  the  computer  has  begun  execution  of  the  code,  the  user 
will  be  prompted  interactively  for  the  names  of  the  files 
containing  the  cross  section  input  data  and  the  problem  input 
data.  The  user  will  enter  the  filenames  at  the  keyboard  and 
the  computer  will  continue  executing  the  program.  Execution  of 
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the  program  will  be  complete  (usually  after  between  1  and  1.5 
hours)  when  the  message  FORTRAN  COMPLETE  is  shown  on  the 
screen  and  the  "$"  prompt  returns. 


D.  SAMPLE  CYLTRAM  INPUT  FILE 

The  input  file  to  execute  an  ITS  run  consists  primarily  of 

order-independent  keywords  followed  by  keyword-specific 

parameter  values.  The  following  sample  input  file  for  a 

CYLTRAN  run  illustrates  the  use  of  the  keywords  and 

parameters. 

ECHOO 
NO  DUMP 
TITli 

391  M«V  Electrons  •  10.0  PMMA.  39cai  Ht.  40.0  PMMA  6  Lf  Arriyt 

. . . 'GEOMETRY . . . . 

GEOMETRY  78 
'Zontl  (PMMA) 

0.00  8.42  0.00  22.0  1 

*  Zmms  2  •  13  (UF  Detector  Arriy) 


8.42 

8.72 

0.00 

0.1 

3 

8.42 

8.72 

0.1 

0.85 

0 

8.42 

8.72 

0.85 

1.05 

3 

8.42 

8.72 

1.05 

1.85 

0 

8.42 

8.72 

1.95 

2.05 

3 

8.42 

8.72 

2.05 

4.85 

0 

8.42 

8.72 

4.85 

5.05 

3 

8.42 

8.72 

5.05 

8.85 

0 

8.42 

8.72 

8.05 

10.05 

3 

8.42 

8.72 

10.05 

18.85 

0 

8.42 

8.72 

18.85 

20.05 

3 

8.42 

8.72 

20.05 

22.0 

0 

*  ZOM  14  mt) 

8.72 

47.72 

0.00 

22.0 

2 

*  Zones  15  - 

28  (If  Detector  Array) 

47.72 

48.02 

0.00 

0.1 

3 

47.72 

48.02 

0.1 

0.85 

0 

47.72 

48.02 

0.85 

1.05 

3 

47.72 

48.02 

1.05 

1.85 

0 

47.72 

48.02 

1.85 

2.05 

8 

47.72 

48.02 

2.05 

4.85 

0 

47.72 

48.02 

4.85 

5.05 

3 
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4772 

48.02 

5.05 

8.95 

0 

4772 

48.02 

8.05 

10.05 

8 

4772 

48.02 

10.05 

19.95 

0 

4772 

48.02 

18.85 

20.05 

3 

4772  48.02  20.05 
*Zmm27  (PMHIU) 

22.0 

0 

48.02 

58.44 

0.00 

22.0 

1 

*  ZMts28- 

38  (Li  Ditictor  Array) 

58.44 

56.74 

0.00 

0.1 

3 

50.44 

56.74 

0.1 

0.95 

0 

58.44 

56.74 

0.85 

1.05 

3 

58.44 

56.74 

1.05 

1.85 

0 

58.44 

56.74 

1.85 

2.05 

8 

56.44 

56.74 

2.05 

4.95 

0 

56.44 

56.74 

4.05 

5.05 

3 

56.44 

56.74 

5.05 

0.95 

0 

58.44 

56.74 

8.05 

10.05 

3 

56.44 

58.74 

10.05 

10.95 

0 

58.44 

56.74 

19.85 

20.05 

3 

56.44  58.74  20.05 
*  Zom40  (PMMA) 

22.0 

0 

58.74 

65.18 

0.00 

22.0 

1 

*  ZMWS  41  • 

52  or  OctectM*  Array) 

65.16 

65.46 

0.00 

0.1 

3 

65.18 

85.48 

0.1 

0.95 

0 

85.18 

65.46 

0.05 

i.(e 

3 

65.16 

85.46 

1.05 

1.05 

0 

65.18 

85.46 

1.95 

2.05 

3 

65.18 

65.46 

2.05 

4.95 

0 

65.18 

85.46 

4.95 

5.05 

3 

65.18 

65.48 

5.05 

8.95 

0 

65.16 

65.46 

9.95 

10.05 

3 

65.18 

85.46 

10.05 

19.95 

0 

6516 

65.48 

19.95 

20.05 

3 

65.18  65.48  20.05 
*  Zmm  53  (PMMA) 

22.0 

0 

65.48 

73.88 

0.00 

22.0 

1 

*  Zqm<  54  •  1 

B5  Off  Oetector  Array) 

73.88 

74.18 

0.00 

0.1 

3 

73.88 

74.18 

0.1 

0.95 

0 

73.88 

74.18 

0.05 

1.05 

8 

73.88 

74.18 

1.05 

1.95 

0 

73.88 

74.18 

1.95 

2.05 

3 

73.88 

74.18 

2.05 

4.95 

0 

73.88 

74.18 

4.85 

5.05 

3 

73.88 

74.18 

5.05 

8.95 

0 

73.88 

74.18 

9.05 

10.05 

8 

73.88 

74.18 

10.05 

10.95 

0 

73.88 

74.18 

19.05 

20.05 

3 

73.88  74.18  20.05 
*  Zmm  66  (PMMA) 

22.0 

0 

69 


74.18 

82.8 

0.00 

22.0 

1 

*  ZoMt  87  - 

78  OF  Oatutar  Array) 

82.8 

82.8 

0.00 

0.1 

3 

82.6 

82.8 

0.1 

0.85 

0 

82.8 

82.8 

0.85 

1.05 

3 

82.6 

82.8 

1.05 

1.85 

0 

82.6 

82.8 

1.85 

2.05 

3 

82.6 

82.8 

2.05 

4.85 

0 

82.6 

82.8 

4.85 

5.05 

3 

82.6 

82.8 

5.05 

8.85 

0 

82.6 

82.8 

8.85 

10.05 

3 

82.6 

82.8 

10.05 

18.85 

0 

82.8 

82.8 

18.85 

20.05 

3 

82.8 

82.8 

20.05 

22.0 

0 

. * . *********80URCC*“ . 

EECTRONS 
BVERGY  381.00 
CUTOFFS  1.803  0.01 
DIRECTION 

. . . OTHER  0PTUIN8** . * . ***** 

RANDOM-NUMBER 

10234 

BiCTRON-ESCAPE 

NBINT10U8ER 

10.0  20.0  30.0  40.0  50.0  70.0  80.0  120.0  150.0  180.0 
HISTORIES  10000 

E.  SAMPLE  CYLTRAN  OUTPUT  FILE 

The  following  sample  CYLTRAN  output  file  was  generated 
using  the  sample  input  file  given  in  section  A. 
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PI06UM  n  -  lEdMia  I,  in? 
VEiaMZ.1 


JMM/LHUim 
TIOMAS  IL  MEKHOM 
IHHUP.  BUSK 

tlUaU  MTHMl  UUMTnSt 


CVITIAM 


***  MUlii  fKCtSM  UM  *** 


MIMKI  V  lUTOlUt 
I 


*  CMIIU  UY  CIOS  8CTUIII  UTA  FM  HUnW  MIMKI  I  * 


MPAI  RTAI  MTU 
Ti  1  AS 


ENFUIFS  (MEV) 


IBOUOMH 

tW.tWflW 

wo.woow 

SAS.Na»7| 

Aouoom 

mnom 

tW.tUQtt 

ISIMIIQIII 

1I0JI000I 

nmm 

5UMm 

4IJI0QN 

muQQi 

ZMSOOM 

lO.IMIH 

II.IWOOI 

Ugggil 

I.IMOW 

^.WHH 

AJOOON 

uum 

unm 

1  jum 

IJHON 

i.wmi 

I.IMOW 

unm 

unm 

umn 

liNON 

LISNN 

L1NNI 

unm 

I.HWW 

I.MlWt 

I.I1WW 

l-tWBW 

unm 

ui5m 

uinn 

mom 

UQIM 

IJMSM 

OHMM 

unm 

uiim 

unm 

TSTAl  ATTEMIATlOil  COEFEKIElin  (CM2/C) 

UITIS 

IJ17I41 

UITUS 

U1717Z 

U1U1Z 

U1I553 

uiim 

M15ni 

IJ15I55 

U1477I 

U14717 

Uiim 

uiaii 

U154a 

IJ1I544 

ui7m 

U2II7I 

tmm 

uziss 

UZIIII 

IJS7II 

imsi 

Uf7f7] 

IJ55II5 

UltSTS 

U7l2t5 

IMIIW 

IWMM 

Limn 

LIIAKI 

L1I22U 

L144S72 

iiilui 

l.1712n 

Limn 

UI774I 

liiua 

liTSAM 

I.5I755I 

IMIWt 

L17niS 

1215552 

15.114517 

52iS4i05 

121215l» 

aiJKiTi 

mmszs  277U15I25 
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UTM  If  SATTgUK  flB  PMI PUHCTIM  Tl  TTAl  ITTUMTIM  CienCEin 


IJ8M0I 

1J0000I 

1J000II 

1JI00I0 

1J00000 

IJOOOOO 

IjgggH 

ijonoi 

9999999 

9999999 

9999999 

9999999 

f  W99W 

0J9SIIS 

1H9037 

1000007 

■  999999 

ojoons 

9999999 

9999992 

9999999 

9999999 

9999999 

IJ99S5I 

1000015 

1000177 

1000700 

im502 

USI737 

1007102 

UOOIOl 

10X0010 

1054010 

1021730 

U545I2 

17I4I4I 

1305241 

1204170 

1IIS454 

1033437 

0J13I20 

1004071 

1001733 

9999999 

1III237 

umn 

lATV  OF  SATIEIIM  TO  SATTEIMG  PIU  PAM  PIOOICTiM  ATIBHATHI  COEFFICIEini 

U14I54 

U1I471 

1024020 

1I3I1I7 

1I3I12I 

1051574 

U7I9II 

1100274 

1101402 

1202014 

1202070 

1300070 

0JI7431 

14557S 

1500550 

1000317 

1700072 

U4IISI 

9999999 

1  999999 

1I243II 

1051200 

1075017 

1003010 

1IIISI7 

K  SEU  HMUTIilll  UTA 

imOINC  EKICY  (MEV),  NOTDEFFECT  EFFtCIENCY  Ml  FIUUEKEIFF  EFFIOEItCY 
I.M0Sa  14l2tn  1.001141 

I  X-MY  ENEIBEI  (MEV) 

LOflOSB  OJOOSB  1000533  1000533 

I  X-UY  ACaMUATEI  KUTIVE  MTENSmEt 
1J00000  1J0000I  1.000000  1JOOOOO 

AOCEI  ELECnOtl  EKICIEI  (MEV) 

IJ00533  0JI0O533  0J00533 

A06EI ELECTHM  ACCOMOUTEI  lElATIVE  WTENST1E3 
1J00000  1J00000  1J00000 


*  GAMMA  lAY  CUB  MCTHM  lATA  Ell  MATEHAl  MMIEI  I  * 


MPAl  RTAI  MTAX 

n  I  41 


HEKEKMEV) 


IMJIQHI 

HUHOH 

jujuin 

m.wwH 

299999999 

isumoH 

101000000 

HjOHH 

HNION 

ll.WWIt 

W  WMM 

21000000 

ISJONH 

11000000 

i.lWHH 

uooooo 

4JHHI 

uooooo 

2JOOOOO 

1JII00I 

1000000 

UONOI 

ISOOOOO 

1400000 

1200000 

1150000 

1100000 

9919999 

t-ftWH 

t-tItllW 

1030000 

1020000 

U15000 

U10000 

uonoi 

IJt4MI 

1003000 

1002000 

1001500 
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iiiilll 


TITM  ATTHMim  nEFFnEKTS  (CMZ/E) 


tmta 

1121405 

M2III2 

uiim 

UllSO 

UIII4I 

U1I454 

U17ISI 

U17127 

UIUK 

UIMII 

Uin4l 

IJ1I331 

U1I454 

U171II 

1III222 

IIZIHI 

1122441 

IJ2S34I 

U27ai 

uniii 

1I3SI14 

1I444II 

1151754 

IJI3SII 

IJ7I7I7 

1111531 

1117171 

UI54M 

IIMSM 

1122757 

III434I 

1I5III1 

1111411 

IITUII 

1115134 

U2I2S4 

UIIKI 

iTinsi 

1J2II21 

1132221 

11211311 

24.711111 

U.7II3II 

11U1I335 

l35J34lt2  1411555122  424U1III5 

MTM  IF  MTTBUK  MS  PU  PIMICIIM  Tl  TTAL  AnBHTM  CeFUEMR 


IJOOOH 

IJOOOH 

1JI00II 

1-WMM 

190IIII 

1911010 

IJONH 

f  f9H9l 

IttWW 

UOUtt 

imon 

4444444 

ttOVIff 

1099917 

UH9W 

1999995 

1009903 

1III1I1 

ItWOff 

1999937 

f 

U999II 

1099972 

4ttow 

I  f9qQ?S 

UOTMt 

1999770 

4949471 

1117151 

1994134 

11SI1I4 

1904021 

U177I7 

1II3I4I 

1759132 

1202291 

1123930 

1037550 

U1903I 

1I07III 

1902321 

1000997 

1000313 

1000140 

1000047 

lATM  IF  SCATTOilU  Tl  SCATTEIMS  PUS  PAI  PUDOCTUR  ATTENUATIM  CIEFFICOT3 

U11744 

U1480S 

1919974 

0924110 

l.tYOMS 

1041343 

1035121 

1131503 

1105703 

1210437 

1257325 

1313035 

■  MMW 

1529450 

1021034 

1743330 

1II334I 

1109751 

1.100001 

1904074 

1037420 

1II750I 

1090730 

1990030 

I  am  BRQATIM  UFA 

IINDIM  HEKT  (MEV),  PHTOEFFECT  ffFIdENCY  Ml  FIIUIESCEMT  ffFICffilCT 
UIOSSS  1.741111  IJ01I4I 


K  l-UY  EKIOEI  (MEV) 

U00S33  IJQ0533  I.U0533  1000533 

I  X-UY  ACOMIUTEI  lEUTIVE  MTEMmES 
1J0000I  IJOOOOO  1J00000  1J00000 

AOOEI EIECFIM  ENEICIES  (MEV) 
tMOVff  IM0533  IM053I 

AOGEI  aECTlOH  ACOMOUTEI  lEUTIVE  MTENSmEI 

ijuou  ijhoh  ijioon 


*  CAMMA  lAY  onss  SniM  lATA  FH  MATEIiAl  MWia  1  * 


MPAI  RTAI  MTAX 
S  1  41 
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BBSKSOIEV) 


IIMIOMII 

nuniN 

mjHm 

SU.I0MI 

^njnm 

3HMM 

mjtMM 

ISUMM 

tNJHON 

tUONH 

lUHiii 

8JHIII 

4tfWWt 

Ttr**** 

lumu 

UMM 

UUM 

uonii 

UHNI 

4JIIM 

UUM 

ZMM 

UUM 

UMM 

UNM 

UUM 

UUM 

L4UM 

UUM 

liHON 

ilSIM 

ilUM 

UUM 

UUM 

usnu 

UUM 

UMH 

U2IIM 

U1SM 

UUM 

UUM 

tfttttt 

UMM 

UUM 

UUM 

UUM 

UUM 

mu  ATTBUIM  CfffnCOTS  (CMZ/C) 

U1M51 

U1K74 

UUM 

U17M 

U17S52 

U17U2 

UINK 

U1I17I 

U1S41S 

U1S114 

U14M 

U14I24 

UI4IN 

U14I2I 

U1SMI 

U1I7I1 

U1M1 

U2US2 

IJ2»S7 

UZS451 

U2tM 

uniu 

U4iia 

U47M 

Losm? 

IJI$44I 

U74S2I 

UUM 

UtSM 

UM41 

1113111 

L1Z453I 

Li4asu 

L1S1I21 

LIUTTI 

LIMU 

liZ7l94 

1.334742 

I.7M744 

1.741117 

uuia 

1UIS241 

vjBms 

12.734351 

213jOM3S 

M12I321  147U344I7  42f1JM21 

uni  If  KJtminc  nn  pu  pudoctim  ti  mu  atterutim  offncHR 

1J000M 

1J000M 

UUM 

1MM 

UUM 

UUM 

i.lMWI 

UUM 

ffgQOM 

1199991 

9999999 

UUM 

f  |<199ff 

UUU4 

9  99WM 

9999999 

UUM 

99999H 

ffaouM 

UMTS 

UUM 

UUM 

US9914 

UUt7l 

UUTU 

UUM 

UMI1 

UUM 

mwff 

US2UI 

U7I452 

US4M 

U94I74 

U2ltl7 

LTinii 

lAUSn 

U1tS24 

imm 

U3HU 

U1SM 

un547 

uum 

UUI5I 

UU275 

UU12S 

UM43 
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UTM  IF  SATTaiK  Tl  KATIINM  PIU  FU  PIWICTin  ATmUlM  OffFICIENIS 


U12147 

U152II 

U2in2 

LNISIS 

1142121 

IJI52M 

UIIISI 

L1S2T 

LimtA 

LZ22I2I 

L2I3M7 

U1IS3I 

LAHTI 

LS35III 

unn 

174107 

umi 

IJ7MI4 

unm 

UMTI 

LUMM 

LNIAT 

f  fffffl 

Lmm 

I  aai  MDITIM  MTA 


IMOM  BBIY  mi).  PHTIEFFECT  EFF KillCY  AM  HHIFtCEIT  TOKY 

iiini?  uzTin  ijm7i 


1 1  lAT  BEKH  (MFV) 

I.I00II7  UUn?  I  WMf7  t  Mflit? 

I X-IAV  ACCMHAia  IBA1NE  HIUaTKt 
I.MOBW  1J000U  l-WflOM  1,MflOM 

AKB  BfCTIH  BEItSI  (MFV) 

unit?  IJIMI7  unit? 

AUFI  BfCTIM  ACCWmUTEI  KUTIVF  HninTlES 
1J900M  1J0000I  tJOOOOl 
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BKIM  cm  SCTIM  MM  FUM  PnCMM  HTMC 


mn  V  sn  M  htvu  MFf  =  i 

. Cnn  tafloni  for  Ml  Heir  Bactnm  la  MIU,  tnaoan 


lUTBUl  won  KTM 
I  iniSK+ll  USIK+N 
Z  A  N 

LIOOOK+II  limiE+ll  U8SUEI1 
unuE+ii  LiaiiE+u  Lam+M 
LUOOK+II  l.1U99E+n  UIIUE+M 

KnmMDPBnnillAlKTCIKKMUX  BUI  EMM  IMAX  UMT 
I  S  I  I  I  1  1  I  M  Uimai+tt  1J27344aE+N  7J13il7K+l1  I 

MTAPIff  MM  FM  UTAFAC  8T  1  M  M  ]  121  ZiltMTE+ll 


AUTEHAi  non  KMn 

2  LlUaf  U  M3SNF+M 
Z  A  II 

I.7I8MI+I1  114II7E-I-I2  UttUE+M 

uum-t-ii  Liasx-t-K  itmie+n 

usTmMDrmnMiKTCMcniiiui  emax  bin  imax  but 

2  i  I  I  X  1  1  I  M  UIHMAE-i-K  1j27M4K+N  iH2l33K+l1  2 

MTANB  MM  FH  MMPAC  ST  2  M  »  I  121  U42t3S+l1 


MATEHAl  won  STMI 
I  imS+11  LnSNE+M 
Z  A  N 

UNME+II  US3SK+I1  liXTSK+N 
UNOS+II  limE+R  I.732SK+N 

MXniMMlZrBERBNilAlKYCICTCMUX  EXUX  EMM  MUX  IMAT 

I  S  I  1  I  1  1  I  M  UliaUS+R  1J27M4S+M  U2I72NE+I1  2 

MTAW  MM  FM  MTATAC  ST  I  M  M  I  121  tJ2]72S+l1 


cailim  /  THAI  S/n  MTBS  FM  MTAPAC  ST  1 
OMiUTIK  MEMBIUMK  cm  SCIWB  FM  MTATAC  ST  I 
CMHATM  HEMCnUSIK  AIKm  IBIWmW  FH  MTMAC  ST  I 
UUKAn  -  EpniMBE  BWMR  FM  MTEIPBATMI 
X  X  Ml  nHKTM  FH  MTAPAC  ST  I 
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purmEcnsN  MKHu  Knuniin 


m  BECTM  BBCY  mflSM  MHIWnM  (lUI) 


*  KEMKMKMPn  * 


78 


OMPUBM  V  mue  KUKMEITII  n  MIKATMt 


mn  9  HAiBuu  H  cm  scnni  ni 

MUNM  MMn  V  EBIEMR  ■  A  nCIlEM  MATEMl 
men  K  BfCnM  CHS  SCIM  BBGY  cn 
aEcnM  BBCY  M  men  » iiwnM  iiEiit  p«nii  SBGT 
nmR  BBCY  oa  men  m  suipum  isMt  pumr  bbcy 
aEcnM  ucu  m  men  m  tuiPUM  BfcniR  ecATmae  ami 
PMm  AMU  eea  men  fab  sampum  eefiit  PHm  amu 
piem  oEieT  eea  men  m  sampum  sbems.  pieni  amu 
BKnM  BBBir  en  men  m  sampum  sbemi  piem  amu 
pieni  Ener  en  men  fh  sampum  peen^aEcnei  iKcnM 
BEcnsi  AMU  en  uMn  fm  sampum  pien  Eucnei  hsecimi 
pieni  BBST  en  men  fh  sampum  pam  oEcnM  OEief 
aEcnsR  BBSY  en  men  fm  sampum  pam  UEcm  bbst 

BABSSUI FWCTIM  Cn  FH  SAMPUM  OECnSN  ENEBSY  ISSS  SnACBUM 

lAMAI  FBMnM  Bn  FH  SAMPUM  OECniN  ENESSY  IBS  SnABBUM 

MAANOM  MMIEB  U  TABUS  9  PienN  CSOSS  SECDOM 

MAXNOM  men  v  picnii  ciiss  SEcm  tabu 

MMBEB  9  BECTHM  ESCAPE  EKICY  SIM 

HMBEB  BE  PIBnM  ESCAPE  EKICY  HM 

MMIEI  9  BECnON  ESCAPE  PIUU  AMU  SMS 

MMIEI  9  mm  ESCAPE  PfllAI  ANCU  IMS 

men  9  tmt  specthm  bbcy  cn 

HMKI  V  PKSE  KKIT  BBCY  HM 
MMIEI  9  BECIIM  BIX  EKICY  IMS 

MMIEI  V  pam  Bn  bbcy  hm 

MMIEI  9  BECnON  BH  OKS 

MMKI  9  PIBTM  BH  OKS 

MMKI  9  PlUmi  OKS 

SOI  9  HniY  MFFEKinAL  HEMS  HSniBIfTIW 

SQE  K  SMCIY  HFFEKimAl  HEMS  HSniBimU 

SDE  9  enBSMtT-SABMEISUI  AMUAI  HSmUTHW 

M.  9  BECTIM  ESCAPE  UWniAI  AMU  HM 

M.  9  PMTni  ESCAPE  AZWCnAl  AMU  HM 


/HUT  = 

I  / 

11 

imat 

=  1 

/  1 

/HUAI  = 

14  / 

H 

/imp  = 

K  / 

K 

1919  = 

41  / 

41 

/UMAX  = 

»  / 

a 

/nAM  = 

21  / 

21 

/HUM  = 

15  / 

25 

/HTAM  = 

1  / 

1 

/■B  = 

11  / 

17 

/MB  = 

11  / 

21 

/MPS  = 

1  / 

17 

/nps  = 

21  / 

21 

/MBAS  = 

INI  / 

1MI 

/MAR  = 

SIM  / 

5HI 

/HTAI  = 

^  1  / 

IS 

/HTAX  = 

41  / 

n 

/UMAX  = 

11  / 

a 

/UPMAX  = 

1  / 

a 

/KMAX  = 

11  / 

a 

/IPMAX  = 

1  / 

a 

/USPEC  = 

1  / 

51 

/USMAX  = 

•  / 

a 

/UFMAX  = 

1  / 

a 

/UFMXP  = 

1  / 

a 

/HU  = 

1  /  1M 

/HUP  = 

1  /  1 

in 

/Mzn  = 

71  / 

1H 

/  ROAM  =  47»  /  4725 

/  MHS  =1115/  UK 
I  MB  =  Z11Z  /  ZtU 
/HMAZ  =1/1 
/KPMAZ  =  I  /  1 
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IHMBEI  K  IHT  2IHS  =  7t 


LffT 

■MT 

MU 

HIB 

BECIlHi 

mm 

PUUE 

PIAM 

UMH 

HUMS 

aiKE 

FMCM 

(CM) 

(CM) 

(CM) 

(CM) 

(MEV) 

Fucim 

1  1  UNMI+II  L42MUE+II  UNME+H  ZiMMI+ll  Unf+N 

{  I  L4ZIMI-t-H  UZflME-fll  UMME-fM  IJMME-II  IJM-t-ll  UNE-t-N 

1  I  L4ZMIE-t-ll  L72IIIE+II  IJNMEII  UMMEI1  IJHE-l-ll  UtK+M 

4  1  ivmt+m  utuK+m  lshmeii  ushk+n  iJis-i-ii  une-i-ii 

5  I  i42aUE-fM  L728ME+H  USME-i-ll  IJaME-l-ll  IJHE-t-ll  UNE+M 
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M  I  IJ44ME+I1  SJ74ME+I1  2MSME+I1  2MMK+I1  IMK+M  IMK+M 

M  1  SJ74ME+I1  U1MK+I1  UMME+M  IMMK+II  IMK+M  UME+M 

41  1  U1MK+I1  U4MK+I1  UMME+M  IMMKI1  IMK+M  UME+M 

42  I  UIMKfll  U4MK+I1  1MMK-I1  UIMK-11  IMK+M  IMK+M 


80 


4S 

3 

ISIIIK+II 

U40002  +  01 

1500002-91 

IJSI002  +  0I 

1J032  +  0I 

19002+90 

44 

1 

LSIIIIE-i-lt 

1.540002  +  01 

1JS0002  +  90 

1J5I002  +  00 

1JI32+90 

10002+10 

45 

3 

LSIIIK+II 

1540002+11 

1J50092+I9 

2.KIIK+II 

1JI32+I9 

19002+11 

« 

1 

L51Nlf+l1 

1540002  +  01 

2JS0002+I9 

4JSI0I2+I0 

1J032+I0 

10002+00 

47 

3 

UlUlf+ll 

1540002  +  01 

4J50002+I9 

US0002  +  0I 

1JI32+0I 

10002+10 

41 

1 

LSIUK-i-ll 

1540002+11 

1050002  +  91 

USMK+OO 

1JI32+0I 

11002+91 

M 

3 

LSiniE+ll 

1540002+11 

1K0002+N 

1JIS002  +  01 

IJi32+9l 

19012+91 

SI 

1 

uinK+ii 

1541002+11 

IjHSIK+II 

IJKOOI+ll 

ijn+H 

19902+19 

SI 

3 

utnK+ii 

1541002  +  01 

1J95I02+I1 

LIIS002  +  01 

1JI32+9I 

10002  +  09 

52 

1 

uimE-i-ii 

1540002+11 

2J05I02+It 

2illl02+l1 

1jn2+9l 

10002  +  09 

S3 

1 

L54INE+I1 

7JSII02+I1 

1900002  +  91 

1211102+11 

IJIS+II 

19002  +  01 

54 

3 

7JtSIIE+l] 

7.41I0I2+I1 

1990002+11 

1JMI0201 

ijn2+9i 

UOH+II 

55 

1 

7JttllE-i-II 

7.4IIII2+I1 

IJMIWII 

lSHOK-11  tJni+H  11102+11 

51 

3 

7JS800E->I1 

7.411012+ II 

1SIII02  01 

IJ5II02+0I 

IM+M 

10002  +  01 

57 

1 

7J33llf+l1 

7.413012+11 

IJSOOK+M 

IJSMOI+OI 

1JI3I+I0 

11002+11 

51 

3 

7J3SIIE+I1 

7.411002  +  01 

IJSOOK+M 

lOSOOK+N 

ijn+oi 

10002+01 

53 

1 

7JSS01f+l1 

7.410002+11 

2JS0002+II 

4JSI002+II 

1JI32+9I 

11012+11 

II 

3 

7J88UE+I1 

7.410002  +  01 

4JS0002+0I 

1050002  +  01 

1J032  +  9I 

10002+01 

11 

1 

7J880K+II 

7.410002  +  01 

1050002  +  00 

IISOOK+N 

1J032+I0 

10002+11 

12 

3 

7.313002+11 

7.410002  +  01 

19SffiI02+00 

tJ05002  +  01 

1J032+II 

10002  +  00 

IS 

1 

7J3000E+I1 

7.411002  +  01 

IJI5002  +  01 

1JISU2+01 

1J032  +  0I 

10002+00 

14 

3 

7J3800E+01 

7.4II002+01 

1J0S002  +  01 

1005002+11 

1J032+I0 

10002  +  00 

15 

1 

7JSSOOE+OI 

7.410002+11 

2J050tf+0I 

1290002+11 

IJ032  +  0I 

10002  +  00 

tt 

1 

7.410002+01 

1210002+01 

1000002+00 

1200002  +  01 

IJI32+0I 

10002  +  00 

17 

3 

L200002  +  0I 

1290002  +  01 

1000002  +  00 

t.nooo2-oi 

1JI32+9I 

10002  +  00 

U 

1 

OJ00002  +  01 

1290002  +  01 

1J00002-I1 

1500002-01  1JI32+0I  UOK+N 

U 

3 

L200002  +  01 

1290002  +  01 

1500002  01 

IJS0002+I0 

1J032  +  0I 

10002  +  00 

71 

1 

L2OOOO2+0I 

1290002+11 

1JS0002  +  00 

1JS0002+N 

1JO32+0I 

10002  +  00 

71 

3 

L2I0002+I1 

1290002  +  01 

1J50002  +  00 

1050002+11 

1J032+II 

10002  +  00 

72 

1 

L200002+I1 

1290002  +  01 

2J50002  +  00 

4JS0002+9I 

1J032  +  0I 

10002  +  00 

73 

3 

L200002+OI 

1290002  +  01 

4J50002  +  0I 

1050002+90 

IJI32+0I 

10002+00 

74 

1 

Ii00002  +  0I 

1290002  +  91 

1050002+10 

1150002  +  01 

IJW+OI 

10002+00 

75 

3 

L2IOOO2+0I 

1290002+11 

1950002+H 

IJ05002+0I 

1JI32+9I 

10002+00 

71 

1 

Ii00002  +  01 

1290002  +  01 

1J0S01K+I1 

1JSS002  +  01 

1JI32+N 

10002+00 

77 

3 

t200002  +  01 

U90002  +  01 

1J9S002+01 

1105002+11 

ijn+N 

10002+00 

71 

1 

0i00002  +  0I 

1.290002  +  01 

2J0S002+0I 

1200002+01 

IJI32+0I 

10002+00 

81 


*  ttIRCf  KMUTM  * 

ana  HtcnoNS 

TK  MAXMM  ana  e«cy  b  ai  jnn  mev 

TK  anu  afCiHR  arifE  non  b  i  jam  mev 

IK  PUm  CniFF  HEKY  B  Uim  MEV 

caouBUTES  K  TK  NHT  avaaa  n  k  tk  ONTEa  k  tk  aEui  (aaio  aam  ue 
1=  ajnaK-t-ncM  t=  ajaaaK+aai  i=  ajaaaa+ncM 

TK  UBBS  K  TK  aEAM  (iBD  aaKE  B  =  BjaaK+n  w 

lEFEiEiia  aBEcnai  Fn  uBaui  aBTHtanaN  b  aEFBEB  n 
TKTA  =  ajm  KBiEEa  PM  =  ajm  kcbee: 

MamaaBEcnaiiu  taota  h  keeieike  bbecthm 

TK  n/uBMn  BOMB  EnMAnt  AK  aaa  n  la  batcki  k  Tm  BBraiEi  eau 
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ElfCTlM-ESAK  BBCT  CUaKATHNI  (MEV) 

&1JIIII  IIUMM  m.71111  m.ll06l  leJMM  IStlllM 
lt7JMII  7121111  KtlM  UnM 


BOIlM-EKiVf  FHJU  jURU  OASanCATMin  (KGIBD 

lUmi  2IJIIII  Himi  1I.MMI  SUMM  TUMI 
MJOM  12UIIII  15UNM  tlUOOM 


aEcnw-EUPE  jUNin  amee  aAsnunon  (kcked 
IN.000M 
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riYsciuiniMis 


afcnM  eamaaii  mi  iuutm  SKity  im  nuesuM 

nna-M  BECIIM  PUNCTIM 

M  cnna  Mume  uttbm  kfucthms 
iKMsnuMR  MB  cuunousiK  x-UY  lOMTA  nuniD 
nfMSHMiMG  Hinnc  moi  v  bibsm  fum  tahutb  mnniM 
PNiM  PMNca  afCTUNS  nuiwa 

lUTBHlM.  1 

aiCTIM  IMK  AT  MAXMni  aUiKE  ENEMY  B  LTSIME+K  (1/011**2) 
I  X-UY  MANTA  MT  FtUmiEI 
ElECnU  NPACT  BNDATMN  NOT  XAMPIEI 

MATEMUNl  2 

EECnON  MIKE  AT  MAXMUM  SOUICE  EKHY  B  il842SE+ 12  (t/aiX**2) 
I  X-UY  MANTA  NOT  FOUOWEI 
BECnON  MPACT  BNOATBN  NOT  UMPIEI 

MATEIIALNO.  I 

EUCnON  lANGE  AT  MAXMOM  MUICE  ENEMY  B  0J0237E-t-02  (G/CM**2) 
I  X-UY  MANTA  NOT  FOUOINEI 
BECTION  MPACT  BNDAidN  NOT  UMPIEI 

ANNMUTBN  MANTA  FUIVIEI 

TI  VOIDME  UTB  B  I.IIOOOOOOOOOOK-t-11 
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at  MM  Bacnm  •  1U  niMil.  acn  Air,  «J  PMMA,  I UF  Ibna 


IK  panfM  B  ia  pocarr  aMPini  iwe  n  fmbi  b  urn  scaot  avbus  tne  pb  uia  b  Aaja 
AVEMS  aacf  BBCV  =  Uiai+K  I  MEV 
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Mn EMK Hucnm  - Ma-MM nmi cbbatq Boma. mhutm hmtm, i lun 

TIMBMISM 

ma  p-n  mHi 

LWK  4  U1M1  t  l.m-ll  1 

1 1  un  mu  MIBBU  M.  t,  1 1 UYI  fum  iutbmi  m.  z.  etc. 

UK+NM  UK+MII  UK+MII 

eiECIM 

mci  p»  UM 

UK-i-MII  UK-t-NM  usn  I 

I  K-Un  EMM  miEMl  M.  1, 1  l-Un  FHM  MATEUU  ■.  2,  ETC. 

LNE-fim  UK+MII  UK+MM 

UTBMIESAPE 

ma  psc  un 

UK  K  4  1.7S  n  I  U4E  K  2 
I  X  un  FIOM  MATEUU  M.  1, 1 1 UH  FIOM  MATEIUI  ML  2,  ETC. 

LNE+NH  UK-t-MM  MK+M  M 


MMKI  AM  EnCT  ESAPf  FMCnUB 


TURMISni 

. BECrm .  PMTM - 

MMKI  EUGY  cum  MMKI  BBCT  CUm 

1.4K+U  1  2.42EI1  I  I4S74  UK-t-H  tUE-ll  I  2S2142 

KflECTMl 

- aECTUR . PMTM - 

MMKI  nEICT  cum  MMKI  HEKT  CUm 

UK+MMUK+Un  I  t.4KII  2111EK  I  I42n 

UTEIUEKAPE 

. aECTUR . PMTM - 

MMKI  HEUY  cum  MMIQ  BBCT  OHn 

S.4KI2  ITJKM  4  M  US+U  ZJKU  1  mTU 
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LMME+MH  LNUE+MH  Z.1I17EN54  Z.1I17ENH 
LMME+MH  LMME+MH  LMME+HM  LMME+M  M 
ZJUK+11  1  -LZIME  I1  ZZ  LUUE+H  Z  U1tK+ll  1 
1UZ1EHZ7  Lias  MS  LNSTENII  U7I7EH14 
LMME+HM  LMME+HM  LMME+HM  LMME+HM 
LUTS  Nil  LNNEMH  ZJ7HEN  I  LNZTEN  4 
LMME+HM  LMME+MH  LHUE+HH  LMME+HM 
LH74EN  I  LNNEMH  1.7434EN  I  L717HN  i 
LNUE+MH  LUOS+HH  LMME+MH  LMME+MH 
4UUEUZI  ZJ7ZZEMH  ilHTENII  LINS  HIS 


87 


UNK+N  umi+ii 
UmE+N  IJtSK+N 
UMK-i-ll  Z.7IS2E+I2 


UtMI+llll  umi-i-lin  UMK+M 
IJ439MM  UNK+MH  7J77SM 
UOBK+NM  UOni+llll  UOIK+li 


ana 


ENEttY  KPSSnn 

(NOBMMQEI  Tl  ME  KKRT  PUIICIE) 


BSKY  KPnmM  (MEV) 


ZOK 

lUTEIUU 

MAS(CM) 

yKIME(CC) 

51 

1 

U117E+M 

LTMK+M 

a 

1 

UMK+M 

7JaiE+l1 

51 

1 

litlK+M 

lilME+M 

54 

1 

t4IMEe 

La4KM 

55 

1 

UMK+M 

M11KI1 

51 

1 

4JMKI1 

1JI5KI1 

57 

1 

UMK+M 

L5447E+M 

51 

1 

U117EI1 

LTMKII 

51 

1 

UMK+M 

iJiaE+ii 

M 

1 

t4IME+M 

La4KI1 

11 

1 

UMK+M 

U27a+I1 

a 

1 

4JMK+M 

1JISK+M 

M 

1 

UMK+M 

t7U2E+e 

M 

1 

U117E+M 

17MK+M 

M 

1 

UMK+M 

7.7211E+I1 

M 

1 

IjtIK+M 

1.2MK+I4 

17 

1 

t411Ke 

La4KH 

H 

1 

UMK+M 

Ml  IK  11 

M 

1 

4JMKI1 

1JI5KI1 

71 

1 

I.MME+M 

t5447E+M 

71 

1 

uanii 

LTMKII 

72 

1 

UMK+M 

1JiaE+l1 

71 

1 

t4114E+M 

La4KI1 

74 

1 

UMK+M 

U27S+I1 

75 

1 

4JMa+M 

1JI5K+M 

71 

1 

UMK+M 

L7iaE+U 

77 

1 

U117E+M 

L7MX+M 

71 

1 

UMK+M 

7.72I1E+I1 

TITU 

TK  ENEMY  CMSnATIM  FUCnW  B  IIOQMSTTKZME+It  I 


PIN  OMI  e-»  Tnu 
UIME+NM  UMK+NM  UIIKHa  UtIKKa 
UIUE+HM  UME+MM  UMK-t-M  M  UMK+MM 
tM7K+ll  I  7.4IMEKI7  UMK^II  1  li7ME+l1  I 
LiaiI  MM  UME  MSt  UNITE  It  17  LIlIK  nil 
UmE+MII  UMK+mi  UMK+MM  UMK+MM 
L4MSK  I  L7MSHM  ZJMSK  I  LM«IE«  i 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
UI7K  I2  7  U578  HM  UMTE  K  I  UMK  K  I 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
USME  HS  UIME  NM  7JIIK  H1I  l.4aK  I2a 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
UiaE  HM  UaK  MM  U41«  Ma  TilME  Ma 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
S.73tKMM  UMK+MM  IJSISMM  Z.IZ7SI4  41 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
IJiai+ll  1  1J37MI1I1  IJHK+II  1  UII7E+I1  1 
IJ44KMSI  UI77EMM  UMZEMS  lltHEHS 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
lilTS  MIl  UtSZE  MIT  UMK  K  7  UMS  K  S 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
tats  12  II  tlTME  Ha  LMITE  M  4  UMTE  M  7 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
IJ47K  a  I  -tiaK  H  K  7J7ta  M  17  t  jllK  M  11 
UOME+MM  UOME+MM  UMK+MM  UMK+MM 
LTMKMa  IJIMEMa  UIHEMM  UMKHa 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
UMK+MM  UMIEMM  t7M2EMS1  TillKMIl 
UMK+MM  UMK+MM  UMK+MM  UMK+MM 
1JMK+e  l  -taiK+M  I  4J7IS+I1  1  1.443K+e  I 
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dUGEimSTM 

(miAiaa  n  me  Hooun  puku) 

ELECTUNS 

ZONE  Moniui  21  zi  ■  n  PUM  nnci  totiu 

1  1  UOOE+M  L42K+II  LIOK+II  2ilQE-(-ll  ZJIOK  H  1E1.177IE  II  (LTOMEIZ  S-1.4MEI1  I 

I  t  L4ZK-I-II  L7ZIE+II  UOflE+M  1JME  II  LIOOIE  M  I7-1.1I0IE  K  I-L40QIE II  1l-t.72IK-l2  7 

1  I  L42IE+II  L7ZIE+II  tJME  II  LSIK  II  Umf+M  M  UmE+M  N  UmE+H  M  UME+M  II 

4  I  ilZK+M  LTZK+II  LSIKII  IJSK+M  UMK-i-ll  M  lillK  n  {7  SJMKM  S4  IJMKII  21 

5  I  t42K+ll  L72K+II  IJSK+N  IJSK+M  UMK+M  M  UMK+M  M  UMK+N  H  UMK+M  M 

I  S  L42K+II  L72K+II  IJSK+M  2JSK+II  UMK+H  M  UME+N  M  4JIIKI4  SB  IJMIM  17 

7  I  L42K+II  L72K+II  2J5K+II  4J5K+II  UMK+M  M  UHK+II  M  UME+M  M  UMK+N  M 

I  I  L42K+II  L72K+II  4JSK+II  SJSK+M  UOIK+il  M  UMK+M  H  UMK+M  M  UMI+H  M 

I  I  L42K+II  L72K+II  SJSK+M  USK-f^M  UMK+M  M  UMK+H  H  UHK+H  M  UMK+H  n 

II  I  L42K+II  L72K+II  USK+M  IJISE+11  UMK+M  M  UMK+M  M  UNK+M  M  UMK-t-H  H 

t1  I  L42K+II  L72K+II  IJISE+II  IJIS-fll  UOIK+II  M  LIIIK-»-ll  H  UNK-f-ll  M  UNK+II M 

12  1  L42K+II  L72K+II  IJKE+ll  2JIS+II  UnK+H  M  IJNK-M  M  UNK-i-M  M  IJMKM  M 

II  I  L42K+II  L72K+II  2JISE+I1  L2IK-^I1  UMK+H  M  UMK+M  N  UMI+II  H  UMK+H  H 

14  2  L72K+II  4.772E+I1  UME+M  2ilK-t-l1  1 JUK  M  M-UOOK  M  17  7J00K-I4  ILI JME  M  H 

15  I  4.772E+I1  4.I02E+I1  UUE+N  1J0K  I1  LIOOK+M  M-IJOOK 14  ILLUME  13  17  L2IIK  IS  IS 

II  I  4.772E+I1  4J02E+I1  1JME  I1  LSIK  II  LOOOK+H  M  UOOK+H  H  LNIK+M  M  LIHK+H  M 

17  I  4.772E+I1  4JI2E+I1  LSIKII  IJSK+M  LMOM+M  ILLUME  H  1L1.1MKU  2LL4IIKM  12 

II  I  4.772E+I1  4J02E+I1  1J5K+M  1J5K+M  LOOHE+M  M  LMOK+M  M  LMOK+M  M  LMME+M  M 
19  1  4.772E+I1  4J02E+I1  1J5K+M  LMK+M  LMME+M  ILLUME  M  51 -LJIOK  M  41-7 JHK  M  H 

21  I  4.772E+I1  4J02E+I1  LISK+M  4JSK+M  LMME+M  M  LMME+M  M  LMME+H  M  LMOK+N  M 

21  I  i772E+l1  4J02E+I1  4J5K+M  5JSK+M  LMME+M  M  LMOK+M  M  LMME+M  M  LMME+M  M 

22  I  4.772E+II  4J02E+II  5JSK+M  L9SK+M  LMME+N  H  LMME+N  M  LNHE+H  H  LMNE+H  H 

23  3  4.772E+I1  4J02E+I1  USK+M  IJISE+II  LMME+M  M  IJMM-M  M  LMOK  M  H  LMME-M  71 

24  I  4.772E+I1  4J02E+II  IJME+II  1J9SE+II  LMHE+M  M  LMOK+N  H  LMOK+N  M  UHK+N  M 

25  3  4.772E+I1  4J02E+I1  1JSSE+I1  LMSE+11  LMOK+N  H  LMOK+N  H  IJHK-M  H  1JHK-I4  M 

21  I  17721+11  4J02E+II  LNSE+11  2iNE+ll  LIOME+M  M  LNHE+N  M  LNHE+N  N  LNNE+N  M 

27  1  4J02E+I1  SJ44E+I1  LNK+M  L2ME+II  2J1NEI2  7-L14HEI2  LL2NKI2  LL52NEI2  S 

21  3  SJ44E+II  SJ74E+I1  LNK+N  1JME  I1  LNHE+N  IL4JNK  M  41 1JMK  N  27 1JNK  H  21 

29  I  5J44E+I1  LI74E+I1  IJHE  II  LSIK  II  LHME+N  H  LNHE+N  N  UHK+M  H  LMNE+M  H 

M  3  LI44E+I1  SJ74E+I1  LSIK  II  1J5K+M  LMOK+N  N5J0HE  M  II  LNHE+N  N5JNK-M  H 

31  I  SJ44E+I1  SJ74E+I1  1JSK+N  1JSK+N  UHK+M  H  LNHE+N  H  LNNE+H  H  LNHE+H  M 

32  3  LHK+II  5J7K+I1  IJSK+M  2JSK+N  LNHE+N  H-7J00KM  H'L71  ME  21  H-7JHKM  H 

33  I  SJ4K+I1  5J7K+I1  2JSK+M  4JSK+H  LHOK+M  H  UBOK+N  H  LHNE+H  H  LNHE+N  H 

34  3  SJ4K+I1  SJ7K+I1  USK+M  SJ5K+H  LNME+N  N  UHK  M  N1JHK  M  N1JI42E  2I  H 

M  I  5J44E+II  5J7K+I1  SJSK+H  LSSK+N  UHK+N  H  LHIK+H  H  LMHE+H  H  UNK+H  H 

M  I  5J4K+I1  SJ74E+I1  USK+H  1.NS+I1  LHOK+H  H  LNNE+H  H  L71ME  21  H  L7IME  21  H 

37  I  SJ4K+II  5J7K+II  1JKE+II  tJ95+ll  UHK+M  H  LMOK+N  H  UHK+N  H  LNHE+N  H 

H  I  SJ44E+I1  SJ7K+I1  1JME+I1  LHS+11  LHME+N  H  LNME+N  HI  JHK-M  H-1JHK 14  H 

N  I  U4K+I1  SJ7K+I1  2JME+I1  L20K+I1  UHK+M  H  LHOK+H  H  UHK+N  H  LNNE+N  H 

41  1  SJ7K+I1  U1K+II  LHK+N  L20K+II  L750K  I2  LL7HK-H  SL4JI0K 12  12  LINK  M  H 

41  3  U1K+I1  LS4K+I1  LHK+N  1J0K  I1  UMK+M  N1J00K  H  ILLMHE  M  NUMK  M  H 

42  I  U1K+I1  L54K+I1  IJHEII  LSIKII  LHME+N  H  LHME+N  H  LMOK+N  H  LMHE+M  H 

43  I  UIK+11  L54K+I1  LSHE  II  USK+H  UHK+M  H  1JMK  M  H-I.IHK-M  42-1 JHK  H  75 

44  I  UIK+11  L54K+I1  1J5K+N  USK+H  UHK+M  H  LNHE+N  H  LNHE+H  M  UHK+M  M 

45  I  UIK+11  LS4K+I1  USK+M  2JSK+N  LHME+N  H  ZJOOK  M  HUNK  M  ILLHNE-M  H 

41  I  UIK+11  L54K+II  2J5K+N  USK+M  LMNE+M  H  LMOK+N  H  UOHE+N  H  UHK+M  H 

47  3  UIK+11  L54K+II  USK+M  SJSK+N  1 JHK  M  M  UHK  M  N2JHK  M  IMJMK  21  M 
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iSIK+ll  LMK+II  USK+M  USK-t-M  UMK+li  M  UME-t^N  M  UMK-t-M  H  UMI+II  M 
li1K+l1  liUE+ll  USK+N  IJKi-i-ll  UWI+H  H  UME+N  »L7IKE-21  « 

LSIK+11  li4K+ll  IJKE-t-ll  UNK+imiJIlll+IIIIIJIMI+NlilJIMI+MH 


BUIS  vomn 

(NomuiQB  n  SE  BOSH  p/umsn 


ZIK 

lUTHttU  71 

0 

SI 

3 

UIB-t-11 

U4S-t^l1 

s 

1 

UIS-t-11 

U4S-hl1 

S3 

1 

U4B-t-l1 

7JIS-t-l1 

S4 

3 

7JIS-hl1 

7.41S-t-l1 

SS 

1 

7jai-t-ii 

7.41S-t-l1 

SI 

3 

7JIS-t-l1 

7.41S-t-l1 

57 

1 

7JHI-t-l1 

7.41B-t-l1 

51 

3 

7J3S-t-l1 

7.41S-t-l1 

SI 

1 

7JSS-t-l1 

7.41S-t-l1 

M 

3 

7JIS-t-l1 

7.41S-t-l1 

11 

1 

7JSS-t-l1 

7.41S-t-l1 

12 

3 

7JME-t-l1 

7.41S-t^l1 

n 

1 

7JSS-t-l1 

7.41S-t-l1 

H 

3 

7JSS-t-l1 

7.41S-t-l1 

IS 

1 

7JSS-t-l1 

7.41S-t-l1 

M 

1 

7.41S-t-l1 

UMI-t-ll 

17 

3 

UME-kII 

L2IS-t-l1 

M 

1 

UME-t-11 

UME-t-11 

H 

3 

lilS-t^ll 

USS-t-11 

71 

1 

L2Mf-t-l1 

UlS-t-ll 

71 

3 

L2IS-t-l1 

UlS-t-ll 

72 

1 

UMf-t-11 

UME-t-11 

73 

3 

L2IS-I-I1 

UME-t-11 

74 

1 

UIS-i-11 

UME-t-11 

7S 

3 

UlS-fll 

UME-t-11 

71 

1 

liNE-t-ll 

UIS-l-11 

77 

3 

L2Mf-t-l1 

UME+11 

71 

1 

L2Mf-t-l1 

UME-t-11 

Bimm 

■  M  PIM  OKI  tsc  ms 

IJIS-t-ll  UKE+II  UMS-t-MHUIIS+NII-lJIHI-MSMJIISHSI 
uaa+ll  L2MI+I1  UMS+imiJMI+IIIIUlHI-t-MIIIJIIS+llll 
ms+ii  Lzis-fii  LznsK  MJusK4M.iiHicauiHiea 
UME+H  iJMSm  IJMI+NNIJIMI+HIMJNS4«l1^mSM11 
1JIS  II  ISMim  UMS+M  M  UNS+N  M  UMS-t-H  M  UMS+M  M 
LSISII  IJHI-t-M  ZJIDSMTUIMIMHiJIHIMHUNSHM 
1J5S+H  UBS-i-a  UMS-fHHUIMI+HNIJIMI-t-NNIJHS+HH 
1J5S+M  LKS-t-M  lJIISM»UIMf  MM-LaaiaiMJMSMM 
LKS+M  Um-l-N  UNS+N  M  UMS+N  M  UNS+M  M  UNS+N  M 
USS+N  U5S-I-M  LMNE-MS1-UMSMaiJIIS-M»1JMSaM 
SJai+M  LBS-i-M  UMS+MHUIMI+IIMIJIIS+MNIJHS+MH 
U5S+M  IJKE-I-II  IJIMI+IIH1JIISM»1JaS-aMIJMI-MM 
IJei+ll  IJttS+MHIJIMI+IIIIIJIIS4NMIJNS+NN 

IJIS-t-ll  LHS-i-ll  UMS-t-M  M  UMS-t-M  »S.ttlS  21  ll-S.42tS  21  N 
2JIS-I-I1  2JIS-fl1  UHS-i-ll  M  UNS-t^N  M  UMS-t-N  M  UMS-hM  M 
UME-t-M  2JMI-hl1  1i24S  I1  2  1  JUS  12  2S^2J7IS  12  a  1  JUS  II  I 
UlS-i^M  1JQSI1  UMS-t-M  IIIJIHE-l-ll»2JIISM»2JINE  MM 
1JME  I1  UME  I1  UMS-t^M  H  UMS-t^M  M  UMS-t^M  M  UIME-kM  M 
UME  I1  USS-t-M  1 JMS  M  M  UMS  M  lEUMS  M  a  1JMS-M  M 
USS-t-M  IJSS-t-M  UMS-t-M  M  UMS-t-M  M  UMS-t-M  M  UMS-t-M  M 
IJSS-t-M  2JSS-t-M  UMS-t^MM2JMSMMUIMEMMUIMIMM 
2JaE-t-M  4JSS-t-M  UMS-t-M  H  UMS-t-M  H  UMS-t-M  M  UMS -I- M  M 
USS-kM  USS-t-M  UMS-t-M  »SJMS  M  II  UMS  M  »2JMS  M  M 
USS-l-M  USS-t-M  UMS-t-M  HUHS-kMM  UMS-t-M  MUHS-fMH 
USS+M  UOS-t-ll  UMS-t-MMUMS-HM»UIMI-MI7UMS-MI7 
IJIS-t-ll  IJIS-t-ll  UMS-t-MMUMS-t-MHUMS-t-MMUMS-t-MM 
UlS-hll  2J0S-t-l1  UMS-t-MMUNS-t-M»2JIMIMM-2JIHEMM 
2J0S+I1  2ilS-t-l1  UMS-hM  M  UMS-t-M  M  UMS-t-M  M  UMS-t-M  M 


mu  tIMSIl  1  1.777SI1  S-L2tM[l1  UJ4MII2  12 
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BKKY  fffCTM  If  muBvina  afcinn 
(MMKi/imr,  muuza  n  m  kkrt  puhcu) 

E(MfV) 

m  JW  -  SI  JM  MS-i-SS 
SI  JOS  I12JSI  7J7E  S  SI 
SIZJIS  -  2717SI  I.IK  M14 
Z717SI  2UJIS  U1ES  i 
ZMJSI  -  ISJIS  1 J3E  S  S 
ISJSI  tSlSS  ZilE  S  I 
15I.SS  -  117JSI  LUI  S  I 
117JIU  7USI  i14I-n  I 
7L2UI  mis  UK  U  1 

mis-  ijm  1J4IK  I 


EHST  VECTU  M  IfflECTB  EECTinB 
(WMKI/MEV,  NSMUDEI  Tl  UE  MCIOEKT  PUTKIE) 

E  (MEV) 

S1J0S  S1J0U  LSE-t-SS 
SIJOS  IILSS  LSE-»-SS 
ULUS  7717IS  LSE-t-SS 
2717IS  mjis  IJK+SS 
mss  ISJSI  LSE-t-SU 
1S.SSI  -1SLSS  LSE-i-SS 
1S.SS  -  117JSI  LSE-i-S  S 
117JSI-  7IJSI  IJK-l-SS 
7L2SI-  mSI  LSE-t-SS 
mis-  1JSI  iJK-t-ss 


HEUT  SECTS  IF  UIESUY  ESSHC  OECTUS 
(SMIEI/MEV,  UIMAUZEI  Tl  SE  MCOERT  PUTTCIE) 

E(MEV) 

SI.US  -S1J0S  LSE-i-SS 
&TJSI-I1LUU  LSE-t^SS 
ULUS  Z7L7UI  LSE-t-SS 
77L70S  -m.SS  IJIE-t-SU 
mjos  IILSOS  LSE-t-S  u 
1SSSI  -1S.SS  IJOE-t-SS 
TSUS-117JSI  LSE-t-SS 
117.su-  7IJUI  LUE-t-US 
TtJIS-  mSI  7J7ESS1 
mis-  1JIS  i.4»n  s 
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MBKAi  KnuHTin  V  TURauTia  Mi  m\ina  Bfcnm 
(MMn/a,  MHUiQB  n  ai  KKirr  piumaf) 


Pa(IEC)=  UM 

1KTA  (lEt)  HUM 


UMI  - 

1UMI 

U4I+M  1 

1UMI- 

MJMI 

ua+M  2 

2UMI 

MJMI 

UK  11  4 

MJMI  - 

MJMI 

UK-11  4 

4UHI  - 

MJMI 

UK  12  S 

MJMI  - 

TUMI 

U4EI2  1 

TUIM- 

MJMI 

U1E-H  14 

tUOM 

IMJIM 

UK+MM 

1ZUIM- 

1MJMI 

UK+MM 

IM.MM- 

IMJIM 

UK+MM 

Mfiui  MnyrniM  a  utbult  hunk  ekthib 
(NMn/a.  MHULQii  n  mi  kseit  punaf) 


ni(K()=  UM 

IKTil  (KC)  IM.NI 


UOM  - 

1UMI 

UK+MM 

lUOH- 

2UOM 

1.7K  H  45 

2UMI- 

MJiMM 

U2E 12  12 

HMM- 

4LMM 

IJ4E 12  7 

4UMI 

M.IOM 

U4E 12  11 

MJIM  - 

TUOH 

LTH-H  1 

7UMI- 

MJIM 

LME  M  IS 

MJOH  • 

12UIM 

US  M  71 

I2UOM- 

1H.MM 

iME+MH 

1MJOH- 

1IUIM 

UK+MM 

HEICY  VECIU  AM  MMUI  KIWmMI  V  aiCTHM  lUMMinB  Ml  eifCTB 
AZNnUl  HIEMM  B  UQMI  Tl  1MJHM  KMEIS 
(MMia/(MErM),  MHUuza  Tl  Ml  numi) 

TKTA=  UM  1UM  MJM  MJM  MJM  MJM  7UM  MJM  IMJM 
I  (MfV)  1UM  MJM  MJM  MJM  MJM  7UM  MJM  IZUH  IMJM 

MUMMIJM  UK+MMUIE+MMUIE+MMUK+MMUK+MMUK+MMUK+MMUK-i-MMUK+MM 
S1JMI1UM  U4EMSI  UK+M  M  UK+M  M  UR+M  M  UK-l^M  M  UK+M  M  UK+M  M  UK+M  M  UK+M  M 
lIUM-ITtTM  l.tS-H  14  UK+M  M  UK+M  M  UK-t-M  M  UK+M  M  UK-t-M  M  UK-t-M  M  UK+M  M  UK+M  M 
ZTlTM-mJM  UK-H  S  t.lTE-M  M  UK+M  M  UK-t-M  M  UK-i-M  M  UK-i-M  M  UK-i-M  M  UK-t-M  M  UK-t-M  M 
mJN-llUM  tJSK  S  UK  Mil  UK-t-M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M 
IMSM-IMAM  2.IK  M  1  US  M  14  UK-t-M  M  UK+M  M  UK -t^M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M  UK-t-M  M 
IM4M1I7JM  Z.7K  I2  I  UK  M  I  UK  M  il  UK+M  M  UK+M  M  UK+M  M  UK+M  M  UK+M  M  UK+M  M 
1I7J»7UM  t4K  K  SL78  R  7  7.tK-M  M  4J7E  M  M  UK+M  M  UK+M  M  UK+M  M  UK+M  M  UK+M  M 
7U»M.1M  4J7f  M  2  7JS  H  1  5.7S  M  I  UK  M  17  UK  M  17  US  M  M  UK+M  M  UK+M  M  UK+M  M 
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II.1PI-1  JN  LIK  K  {  U7f  K  {  UK  n  4  1421 N  4  IJKIl  S  4JK  M  I  UK  K  14  UK+H  M  UK-t-M  H 
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EKKY  fff nU  AM  AMItAI  HmHIMB  ■  BICIMM  lUnima  AM  UlECTB 
A2NITIAI  HIBNAl  B  UMH  11 IHJNN  KttEES 

(MMKi/ounrB),  MMMAuza  n  m  paihoe) 


TKTA=15I.NI 
E  (MEV)  1HJM 

aiJMI  SIJNI  IME+NM 
SIMM  IIUMI  UK+HM 
312MN  27L7MI  UK+H  H 
27I7MI  mJMI  UK+HN 
ZS4JNI  IKMH  IME+MH 
ll5MH  -ia.4MI  IME+HM 
ISLAOM  •  117JMI  IJK+NM 
117JMI  TUNH  UK+MH 
TUMI  MINI  LNE+NH 

aim  ijm  ime+nm 


ENEMY  VECTU  AM  ANMUl  HSnMtfllMS  V  aECIUH  UIBAUY  EKANK 
ATWniAl  HIEHAL  B  UOm  11  lUJNM  KOaS 
(IIIIMIEI/(MEV*B),  MMALDEI 11  ME  YANCIE) 

iKiA=  im  iMM  2im  am  am  am  tun  am  lam 
E  (MEV)  am  am  am  am  am  tun  am  lam  lam 

Btm  aim  lue+m  m  lne+m  m  imi+m  m  lme-i-h  m  ime+m  m  uk+m  m  mk+n  m  UK-t-a  a  UK+a  a 
aimiitm  LaE+aauiE-t-aauaE-t-aaiJiE+aauK+aaiME-t-aaiJK+aaiME+aaiaE+aa 
iiLm37L7a  uiE+a  a  ua-t-a  a  ua+a  a  ua-»^a  a  ua+a  a  ua+a  a  ua+a  a  ua+a  a  ua+a  a 
m.7Mmm  ua-t-aaiaE-t-aaija-t-aaija-i-aaiJa+aaija+aaua-i-aaua-t'aaija-t-aa 
mmiam  ica+aaua+aaua+aaua-i-aaua+aaua+aaua+aaua+aaua-i-aa 
iamia.m  ua+aaua+aaua+aaua+aaua+aaija+aaija+aaMa+aaua-i-aM 
ia.mii7m  ua+aaua+aaua+aaua+aaua+aaua+aaua+aaua+aaua+aa 
ii7JM7iia  ua-t-aaua+aaua+aaua-i-aaua-t-aaua+aaua+aaua+aaua+aa 
712M  aia  ua+a  a  ua+a  a  ua+a  a  ua+a  a  uaa  a  ija  a  a  i.ia  a  a  ua+a  a  ua+a  a 
a.iM  im  laE+a  a  4.7a  a  «  in  a  iz  liie  a  7  ua  a  11 1.7a  a  1  zn  a  u  zm  a  71  ua+a  a 
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non  VECTU  «■  umu  MnumaB  if  Bicnw  uibhuy  e»pm 
UIMTUI  ITBViU  B  UNM  Tl  lOJIlll  KCKES 

(MiMia/(MEira),  mmuiqb  n  me  rumn 


TKTA=1SIJM 
E  (MEV)  HUN 

N1JIN  K1JIM  UK+NN 
SUM  - HUM  UK-i-HH 
IILHN  ZTITM  UK-l-HH 
ZTlTM  -aUM  UK-t-HH 
ZMMI  ISMI  UK-t-NH 
tSMI  •  ISI.NN  UK+N  H 
1SUM  -1I7MI  UK-l-NH 
ll7Ma-  7IMI  UK-t-NH 
TUM-  MEM  UK-i-HH 
H.1M-  IMI  UK-i-NM 
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APPENDIX  B.  INTERACTIVE  DATA  LANGUAGE  (IDL)  DATA 
MANIPULATION  AND  PLOTTING  PROGRAM 


IDL  is  a  programming  language  commercially  available 
from  Research  Systems  Incorporated  that  is  capable  of  being 
used  in  a  wide  variety  of  computing  and  graphics 
applications.  In  this  work  the  IDL  language  was  used  to 
write  a  relatively  user-friendly  program  that  performs  basic 
functions  necessary  to  extract  data  from  the  standard 
ITS/CYLTRAN  formatted  output  files,  rearrange  the  data  into 
proper  data  structures  for  plotting,  and  plot  the  data  on  a 
selected  device.  There  is  no  substitute  for  reading  the  IDL 
user  manuals  as  an  education  in  how  to  use  the  language, 
however,  the  following  is  a  user's  guide  to  using  the 
analysis  and  plotting  program  written  for  this  work. 

A.  RUNNING  IDL  AND  EXECUTING  PROGRAMS 

The  IDL  compiler  is  an  interactive  program  that  is 
installed  on  the  VAXStation  3200  in  the  directory  IDL. DIR. 
All  IDL  commands  are  compiled  as  they  are  entered,  whether 
they  are  entered  directly  into  the  compiler  from  the 
terminal  or  as  a  batch  stream  from  a  command  file,  or 
program.  In  this  respect,  then,  IDL  runs  in  a  mode  very 
similar  to  PC  Basic.  To  invoke  the  IDL  compiler  after 
logging  on  to  a  user  account,  type  "IDL".  Some  licensing 
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information  will  be  printed  on  the  screen  and  the  prompt 
"IDL>"  will  appear.  At  this  point,  IDL  commands  may  be 
entered  and  compiled  directly  from  the  terminal.  The  point 
of  entering  commands  this  way  is  that  the  powerful 
capabilities  of  the  language  (which  can  be  appreciated  by 
reviewing  the  commands  and  what  they  can  do)  can  be  used 
without  a  lot  of  programming  effort. 

A  complicated  series  of  IDL  commands,  such  as  the  ones 
necessary  to  rianipulate  the  large  CYLTRAN  output  files, 
should  be  typed  into  a  command  file.  The  format  of  the 
commands  is  the  same  as  if  they  were  entered  interactively, 
but  they  are  put  in  a  text  file  using  the  text  editor 
available  with  the  VAXStation  3200  (called  EVE) .  After 
constructing  the  command  file,  save  it  with  a  name  ending  in 
the  file  extension  ''.PRO'',  for  example  MY_PROGRAM . PRO . 

To  run  an  IDL  program  saved  in  a  text  file,  type  the  IDL 
executive  command  ".run  MY_PROGRAM. PRO" .  The  IDL  compiler 
will  compile  and  execute  the  program. 

Output  from  the  program  can  be  displayed  on  one  of 
several  devices  with  great  ease.  To  display  a  plot  on  the 
monitor,  the  command  "SET_PLOT,  'X'  "  s  ould  be  entered 
prior  to  the  PLOT  command,  which  generates  the  plot  itself. 

X  (for  X-windows)  is  the  default  output  device.  To  plot  a 
drawing  on  a  Hewlett-Packard  plotter  which  uses  the  HPGL 
language,  enter  "  SET_PLOT,  'HP'  "  prior  to  the  PLOT 
command . 
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B.  RUNNING  PLOT. PRO  TO  ANALYZE  CYLTRAN  OUTPUT 

The  name  of  the  IDL  program  written  to  analyze  the 
formatted  output  generated  by  CYLTRAN  in  this  work  is 
PLOT. PRO.  The  main  program  calls  the  three  procedures 
(which  is  the  IDL  name  for  a  subroutine)  called  READ. PRO, 
PLOTDATA . PRO ,  and  REFDOSE.PRO.  The  IDL  command  code  for 
these  four  command  procedures  are  located  each  in  separate 
text  files  of  the  same  names.  To  run  the  program,  all  four 
of  these  files  must  be  in  the  user’s  directory  (they  don’t 
actually  have  to  be  in  the  user’s  directory,  per  se,  but  the 
user  must  be  able  to  tell  the  IDL  compiler  where  to  look  for 
them,  and  the  IDL  compiler  must  be  able  to  access  them) . 

To  execute  the  program  PLOT. PRO,  at  the  IDL>  prompt, 
type  the  command  "  .run  PLOT. PRO  ”.  This  assumes  that  all 
four  of  the  program  files  are  in  the  user's  root  directory. 
IDL  will  compile  and  execute  the  plotting  program. 

As  the  program  is  executing,  it  will  prompt  the  user, 
either  through  the  terminal  command  window  or  by  a  mouse 
menu  window,  for  input  or  choices.  If  the  user  knows  the 
name  of  the  CYLTRAN  output  file,  the  answers  to  some 
elementary  questions  about  the  problem  geometry  and  CYLTRAN 
input  zones,  what  kind  of  output  he  wants,  and  what  device 
he  wants  it  on,  the  use  of  the  program  should  be  self- 
explanatory. 

The  data  reading  portion  of  the  program  is  generally 
applicable.  That  is,  the  PLOT. PRO  procedure  should  be  able 
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to  read  most  of  the  data  usually  of  interest  from  any 
CYLTRAN  output  file,  regardless  of  the  problem  material, 
geometry,  of  number  of  input  zones.  This  capability  is  a 
tremendous  asset  in  itself,  because  the  user  has  the  ability 
to  halt  the  program  after  it  reads  the  data  and  to  enter 
commands  interactively  to  plot  or  manipulate  the  data  in  any 
manner  desired.  However,  if  the  program  is  used  for  graphics 
output  on  any  CYLTRAN  problem  that  deviates  even  in  a  small 
way  from  the  geometry  of  this  work  (cylindrical  with  annular 
arrays  of  TLD's  radially  dispersed  at  regular  longitudinal 
intervals)  the  results  cannot  be  guaranteed.  Proceed  with 
caution. 

More  comments  on  the  specific  use  and  purpose  of  each 
portion  of  the  program  are  included  in  the  source  code 
documentation  for  the  program  itself.  The  user  is  encouraged 
to  read  through  the  source  code  briefly  before  running  the 


program. 


COMMON  PLOT, 

■•duln.iNn,Piiin,ilactaiiM.nrin,s«4Ml  ,s«riuZ,s«riuS,S 
«ariN4,sepiM5,strlM8,ttrw2.ilNt_MNmr.tar|it,Ptf,Kl1.i^lZ,l 
ptl3,rtM,KMts,acctMx1  .acctptyl ,  ii6c«ptx2,acMpty2,$ 
>cci(po3,»cc»pty8,«cctptt4,ioeMrt»4,ntirt,nic>ilci,clnlct,l 
ptotcMM,paiMHi,BixvaLBlniLni>caldNiM2,ttnKl,S 

■llH.Lik.ZMM 
iiMt=0  itlUi=" 
dNte=0 

■tml = {nic,iMi*ir:0,Hat:0,zton:O.O.iPiMK:0.0,rti:O.Q,pa«t:0.0.$ 
■lMrt:0.0.plMtcM:0.0,Ma«s:0.0.wkMM:0.0,MHirBMT(8)} 
strae2  s  {aiuMM,MHlMP:O,Hit:0.ikR;0.O.zri|kt:0.a,pk:O.O.r«it:O.O,S 
•iart:0.0,piNtcit:0.0.Han:0.0.vslHM:0.0,«Mriy;fltarr(8).  I 
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n*is:0.IU:0.0.ilact:0.0,iui:mirr(8)} 

;Tkii  b  dM  nil  mi  sKtiN 
lta|=0 

wHi  lib  EQ  0  «•  iHb 

ciNic«=wmi(rin  CMPMCS  ROUTM'.’Pnccss  a  Ml  H  iataM 
Timca  a  ataiiart  piat','Exll1,TUa=0) 


2:  phtiati 
8:  IM^I 


priM,7ypa  mtna  ta  iai  tka  praraa ' 
prM.'  ar  aitar  btaraetiw  IN.  camaadt/ 
atap 
aid 

PROraaddata 

;Tiia  pracadira  akoiH  ka  abia  ta  raid  aiy  aatpit  Oa  pradicad  by  tka 
;CVITIUN  rre  cada.  k  ny  avia  ba  abia  ta  rtad  JUXOT  aitpk  Baa. 
;bit  ka  aMHy  ta  da  m  Ina  lat  baai  taatad. 

aMMONPtOT, 

■*dnlBx.iaBaa.ndhtt.zfiet  taHit2.aaHatS.I 

aariaa4,aariaa5,aariaa8,atrac2,alact_aiany,tarfat,ra(,i^,pal2,S 

pal3,ral4.ralpta.accaptx1  .accaptyl  ,accaktx2,accakty2,t 

acta|tx8,acca|ty8.accaitt4,accaity4,ntart,8adcbalca,cbalca,t 

iiatchaiea,ranai.nnalainilllNCMl,cfeaica2.atrie1,| 

Iai,tkla.l8iian.cird,tan»tric.iaaatan.tant,biwnk.tan2,8 

allai,UNnM 

■HnNs"  ttkiBs"  tzaaaa^O  ii=0  U=0  ik=:0  tallai=0 
caN=” 

prbrt.  Typa  tha  am  af  tha  m  aatpk  Ih  ' 

priat.  Iamt=  '('  a.|.  Qaiaa.dbntary]cye0.ait;1  M) ' 

rniWumm 

apaip,  1,  Baian 

far  1=0,14  da  raadf,1,card 

WHNE  NOT  E0F(1)  DO  BEGH 

raadi,  1,  caN  ;ba|bi  rndbii  data 
jInIi  top  toy  tosA^y 
I  atnaid(capd,O.S)  ap  '1 . On  bapia 


.Had  tha  aiMbar  al  laaaa  aid  tba  zaaa  piintry  data 
nl  Kd4  dm 

I  atpaid(capd,5,24)  ap  'GEOMETRY  DB>ENDBirT  HPUT  I 

Ifeaibaila 

raadi,  1,  card 
^raad  tha  ■aibar  af  zaaaa 
raadt  1,  lanBat='(24i,l5)',  ma 
nia = raplcata(atrac1  ,iaBM) 

;craita  a  atricdra  far  each  bnt  nia 
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tMpMTKs  {nMtMp,iHta>:0.0,Mt:0.0,iWLlMI,f 
irf|M:0.0,i«:OArMa:0.0,>icat:0.Q,ilNtatt:0.O  ) 
nwtNw  B  Mplnti(twi»ttric,aMwt) 

;1H>  ttcUN  narit  k  zmm  iNMtpy  data  iar  aiy  Mafeap  af  zaiaa. 
;Tlw  Baabap  af  naaa  k  paad  kaai  tba  ITS  aatpai  Ba. 
kp  1=0,2  ia  paaA,1,capi 
paadf,1.zaaataaa 


iaaa.ikfi = naataav-ikli 
laaa  jpkM = zaaataaa-tPiibt 

iaaa,pk = naataav-Pk 


iaaa.alcat = zaaatBaai.alcat 
iaaa.alMtcat=iaaataaa.alMtcat 

;IM  tta  tida 

I  atpadl(caPd,0,S)  aa  '  *  '  tbaa  kaik 
paadf,1,capd 

I  ttradd(capd,0,S)  aa  '  *  '  tbu  $ 
titk= stparid(capa,S.74) 


;lki  tha  aaapfy  data  aad  ptad  I  kta  tia  atpactapa 

I  (atPtPki(straiid(capd,1.B8),1)  aa  'BliRGir  DEPOSmom  $ 

aad  (allai  EQ  0)  tbaa  baik 

;Tbk  aactka  ptada  k  aaeny  depaaltka  data  kP  800-848  naaa. 

;Ta  Ptad  kp  napt  ap  kwtp  naaa,  yaa  aant  cbaaia  tba  ain  af 
,lba  atpactapa  accapdkify.  La.  ta  paad  858  naaa,  cat  tba  wka 
;af  'kwaadt'  ta  850  ap  ,k  laaapaL  R  abaak  ba  r,  tba  aaapaat 
;kwap  aadtkk  af  50.  Msa,  tba  appap  kai  aa  J  k  aataandcair 
;aat  apaal  ta  that  taakap  dMdad  by  50  akan  1,  b.|.  850/50-1, 

:aa  1=0,12 

taawatpac1  =  {na8teak1,aaarikp:0.0,ant;0.0,an8a:0.0,«Bbaaa:0.0,  S 
aaap|y:IRapp(8)} 
taani  =papflcatB(teiapatpac1,S0) 

ppkt, . . * . * . . 

ppkt. '  Iba  papaawtap  'kwaadt'  abaald  ba  aat  ta  tba ' 
ppkt, '  aaaraat  kwap  aadtkk  af  SO  ta  tba  aaariwp  af ' 
ppkt, '  kpat  naaa  k  tba  prabkw. ' 
ppkt, ' ' 

ppkt,  kPBnt=  '('  Typa  tba  vaka  af  kwaadt:',!)' 
paad,  kwaadt 

taak2 = paplcatB(tBBWctpac1,naaa-kwaadl) 
kP  1=0,2  da  paadf,1,capd 
kp  J=0,lk(kwaadt/50.0-1.0)  da  bapk 
paadf,1,taBMi1 

naaQ'50;|'50-i-48).aBakap=taaklM*ap 
naaa*S0r50-i-48)JMt=taaklJBat 
naaQ*S0r50-»-48)  jnaa=taak1  naaa 
nBaa*S0;l*50-»-40).vabaBa=taak1.VBbaBa 
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mtO*50j*50+49).wwfBf=Hf1.Mwriy 
tor  1=0,8  it  ptai(,1,cari 

Niltr 

rtai(,1,tM|iZ 

itM(to«nM*:ZMM-1 )  Mtotr = tMp2.MBitr 

itit(h«nMt:itMt-1)  JMt = tMp2.Mt 
zNt(toiNMt:ztMS-1 )  jnst = tta|i2  jnss 
ZNt(lMMMR:itMS-1).MtoM=tllip2.MlM 

HM  ( WinMH.ZtMt’ 1 )  .NiPijf  =  li■P£«MtP|J 
tlta|=1 


BIDWWLE 

irkt  'M  i(  itta  III  mcher 

rttara 


PRO  ptotiata 

9 

;T1iB  artctiart  It  larinta  ft  prtiact  ptots  tor  tfea  tptdk 
;CyiTlUN  laaaMtry  httap  txwtai  to  tha  Otsto  by  RN  Yaw.  R  wl 

rPTvmqf  Ni  prwnci  MiiMiiini  iripHci  up  Wmjf  iuNr 

;ctaltourattoa,  atbaapb  aay  itbar  prapraai  tor  ptotttop 
;fTf  aatpat  wl  pnbabiy  ba  tbalv  ta  tbh. 

QMIMONPLOT 

■aiatos.iMMs.ndbtt.zfacLafe'iii.iirtot.tartatl  ■artoi2.iwtotS.8 

AWWft  ^ai^ifwa  awa^wa  awA^wa  awaa^V 

aartoal,aartoa5,8artoa8,atrBc2,atoct_aatriy,tartat,ral,raf1,ral2,l 

ral3,raM,ralpt8,accapni,accapty1,accapti2,accap^$ 

accaptK3,accapty8,accapti4,accapty4.iatart,pacbaica,cbalca,8 

ptotrbatoa,raianBi,wuwtwtovatlit,caal,cbalca2,8trae1,8 

flaiittto,lltaaaw,cari,tawpatrac,iaaatawp,lawp1  ,towwaR,tawp2,l 

alto|,Ub,zaaa 

;Tba  parawatar  waiuka  awat  ba  aat  ta  iaacriba  tba  wwtoar  af  aata 
tol  aaaalar  rtopa  aaai  to  tba  prabtow,  La.  8  (78  zaaaa), 

;88  (750  zaaaa),  80  (828  zaaaa),  55  (860  zaaaa),  ate 
ii=0  lii=0  tcbaica=0  Rcbaka1=0 

PPtot, . . 

prtot,  torwat=  '('  Typa  tba  wwtoar  al  Mactar  S 


Z=fltarr  (zaaaa) 

ralaa=0.0 

ztoct=0.0 

ato|ap=0.0 

aartoa=rap8cata(atrac2,wi*iii) 

aartoal  =  aartoa  Iaartoa2=aartoa  Iaartoa8=aartoa  Iaartoa4=aartoa 
aartoa5saartoa  ftaartoa8= aartoa 

prtot, . . . . 

prtot  torwat=  '('  Wbat  to  tba  avaraia  atoctraa  aaany  bl 
MaVTM)' 
raaiL  atoct  aaarpy 
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pr«.' 


;TWi  stctiM  takM  mu  •!  Um  tariet  ■attrtab  bitai  Mti  Mi 
;MsigM  raiMt  tt  dM  iMM  tfetikMi  by  wptrkMit  at  BatM-MIT 
;m  that  a  rtlaraaca  Im  cm  ba  ylattai  aiaiast  ITS  calcalatiww 

0 

tanct=0 


istart=0.0 

raf_iaM 

;***** . . . * . . . . . laSawbu  bap  aaaipa  um  iau  ta  8 

aaparata  atPMUm; 

;aM  far  aacb  pa*a  at  ntbicb  Icm  b  auasapti.  Tba  Maa  I  caMta 

.•tbraapb  tba  Mabar  af  bpM  znm  praaaat  at  aacb  paAu.  Ibb  b 

;«n8Pait  br  aacb Mawtry.  Tba  appar  bat  all b 

;(iaaM-Maiaitarai  zaaM)/#aaaabr  rbpa  •  1 

;Far  aarba  af  8  LT  iatactara  witb  m  btanaaAata  tarpat  laaM  1=0,5 

;Fap  758  bpat  zaaaa,  m  LT  iatMtara,  1=0,62 

;Far  828  bpat  naaa,  aritb  8  LT  iatactara,  1=0,68 

;Far  880  bpi  zaaaa,  wtb  5  Lf  iatactara,  1=0,54 

Ma  aay  caaa,  tba  appar  baR  m  i  b  auiaba-l 


|cacbaica=* 

prbt,  . . * . . . 

prbt '  lira  tbara  aay  tarpat  aectbM  bativaN  iatactart 
rbpa  far  wbicb' 

prbt,  fanMt=  '('  tba  iau  wb  Mt  ba  pbttair  y  ar  a:*,$) ' 

raaC  paacbaicc 

prbt,' . * . **' 

far  i=0.aMialea-1  ia  bapb 
tar  J=0,N_TMS(zaM)-1  ia  bapb 
8  (paacMca  ap  T)  tbaa  bapb 

aarbs1(D.0=zaM(1  -»-n3).0 
aarba2a).Q) =zaM(3-)-i*13).Q) 
acrba3a).® = zaM(5 +1*1 3).0 
aarbs4(D.(D=nM(7+l*13).(D 
aarba5a).a)=uaa(0+n3).(D 
aarbs8(D.(D=zaM(11  +I*13).(D 
aarba1(Di=Z(1+i*13) 

8arba2(Di=Z(3+i*13) 
aarbs3a)i=Z(5+i*13) 
aarba4a)i=Z(7+n3) 
aarbs5(DJ=Z(0+n3) 
aarba8(DJ=Z(11+i*13) 
aaN  aba  bapb 
aarbal  (D.® =zaMO*12).0 
aarba2®.®=zaM  (2+1*12).® 

8arbs3®.® = zaM(4+ri2).® 
aarba4®.®=zaM(8+l*12).® 
aarbs5®.®=zaM(8+i*12).® 
aarbaS® .® = zaM(1 0 + 1*  1 2).® 
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tarin2(DJ=Z(2+nZ) 

sfriMSa)i=Z(4+n2) 

s«rin4a)i=Z(8-fn2) 

Mrin5(DJ=Z(8-«-n2) 

striM8a)i=z(io-i-n2) 

nMm 

nrin1(D.zlact=(MriKia).MU/tirinia).wlHM)*  I 
(ssriMi  (D  JriiM-strinI  (D  JMt) 
MriMl(D.8m=8triw1(D.wtny/Mrt»«1  (D  Jun 
striM2(D.ifact=(MriM2(D.Batt/SiriM2(D.wiM)*  $ 
(strlu2(DJrm-Mrlu2(DJM) 
»tp|ii2(D.8m= 8triit2(D.Mnrp/nriii2(D  JMM 
MriM3(D.ztact=(uriMSa).nss/MriM3(D.MkM)*  S 
(sarln3(DJri|tt-strin3(D.iMt) 
MriM3(D.i8M=MriM3(D.MNny/t8rinSa)Jntt 
s8riM4(D.zlact=(s8rin4a).atM/s8riN4a).wlBM)*  $ 
(tiriM4(D.iri8M-t8rlMl(D.ililt) 
8WlM4(D.<Mt=t8riMl(D.mpp/MrlM<l(DJMti 
tM‘in5(D.zlact=(scriBs5a).MM«/ttPluSa).MlBflw)*  $ 
(tepies5(D.zpltfrt-urksS(0.iiin) 
MPies5(D.<«M=Mri»t5(D.Mtrp/MflM5fl).m 
siriB88(0.ztact=(seriBs8(D.MSt/scpiu8(D.volHM)*  $ 
(s«ries8(D.zPi|ht-siriis8(D.iiin) 
t«'iM8(D  .dou = seriet8(D  .•Mny/Mries8(D  .mss 
sidfop 
sidisp 

sspMI  (0)  J = sspiesi  (0).zlsct 
stpiu2(0)i  s  septes2(0).ifact 
sspiis8(0)  J = ssples3(0).zlsct 
sspiM4(0)i= scpies4(0).zfsct 
tiplss5(0)  J = sepies5(0).rtsct 
sspiss8(0)i=sepies8(0).zfsct 
top  l=1,Mdttles-1  dt  hcdii 
stphsl  (Di = sspksl  (i-l )  J+sepiBs1(D.zfact 
sspiss2(0i=sepies2(l-1)i+stPtos2(D.ztoct 
ssphsSCOi = sspks3(l-1  )i + sspiBsS(D.ztoct 
sspiss4(Di = scpiBs4(i-1)i  4- ssptos4(D.ztoct 
ssPtos5(i)i=sspissS(i-1)J4-ssptosS(D.itoct 
ssptos8(Di= sspks8a-1  )i  ssptos8(D.ztoct 
tidfsp 

ssptos1.psdins=0.0 
ssptos2j>sdhs=1.6 
saptos3j‘tdlM=2.0 
ssptos1jsdtos=5.0 
ssplts5Jidlus=10.B 
SSPiBS8j>SdM=20.0 
hsp:  tonics  sMflWHd'Sskct  a  ptot 

70181  doss  -  al  psdl  Z  to  i/cnZ*,  $ 

Total  doss  •  y  pad(  Z  to  ooM 
'Bpsakdom  by  tooctpoi  typo  -  aao  padbis',$ 

'Stop  ay  Maks  a  Caston  Ptor,'Exn.  TMo=0) 
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;  nfl  wM  Z  ihN  k  riM  ptr  Minrt  mtkwttr 
;  aliw  a  cbaica  al  aatpit  rapMci  iavicaa 
pktdwlca=iwaaa(| 

I'OMaaa  a  mphica  dtvica  -'.XaMawa',  I 
'HPa'.'Pa',Taatacrkt'].TUa=0) 
caaa  pktchaica  af 
1:aat_|lat,T 
Z:  aat.pkt.'HP' 

8:  aat_pkt,'Pa' 

4:  aat_pkt,'PS' 
aadcaaa 

;  pkt  tia  Bataa  4ata  far  rafaraaca 

I  (tariat  aa  5)  thaa  bagk 
.ikara  ara  aafy  8  rafaraaca  fewa 
;  livea  k  iata  far  thk  tanet 
pkt  k,r8f2.z,rcf2.daac,titk=8trtrki(tltk),S 
ytitk='Dos8  (MeV/li|)',$ 
xtkkslaaiptiidkal  Dktaaca  (i/caOD'.S 
]iraB|c=[1.08-08,1.0a-i-02],  t 
araapa = (0,aMx(rcf2.z+ 1 0)],pa]fai = B 
apkt,rcf2.z,rcf2.4aaa,cakr = 8 
xyaata,rcf2.i(rcfpt8-1 ) + 1 .5,1 

raf2.4aaa(rafpta-1  ),'r = Z.Ocai' 
aadf  aka  bepk 

pkt>,raf1  .z,raf1 .48a8.tna = atrtrki(tlk),$ 
ytIks’Dasa  (MaV/kDM 
idUaslaapItadBat  Dktaaca  (|/cai2)',$ 
yraa8a=[.0000001,aMX(r8t1.48sa)],  I 
xraagas  [0,anx(r8f1  j[4-10)],payBi=8 
apktrtfl  .z,rat1 .4aia,cokr=8 
iy8Bts,r8f1.  >afptt-1)-«-1.S,8 
rct1.i8sa(r8fpts-1),'r=0.0cai' 
8pkt.r8f2.i,rttZ.4aa8,payai= 6 
apkt,r8t2.i,ratZ.iaM,cakr = 8 
I  (tarpat  ap  4)  tfeaa  I 
iyaats,raf2.z(rafpta-1 )  -f  1 .5,1 

ratZ.4aaa(ratpta-1).'r=2.0cai'  $ 

aka  8 

iy8ats,rat1  .z(ratpta-l) +1 .5,1 

raff  .daaa(rafpta-1  ),'r = 1  .Gcal' 

aapaka 

apkt,rat3.^rat3.4aaa,p8yai = 8 
8pkt,rat3.z,rat3.4aaa,cakr=8 
xyaata,rat8.z(rafpta-1 )  -f  1 .5,8 
rat8.4a8a(rafpt8D,'rs5.0cai' 
apkt,rat4  j,raf4.iaaa,pa]iBi = 8 
apkt,raf4.i,raf4.iaaa,cakr=8 
aikata,raf4.z(ratpt8-1 ) + 1 .5,8 


rtl4.ioM(rtlpts-1  )/r = 20.0ai' 
x= [10,18] 
y=[8.0«-08,6.0t-08] 

I  (tanit  •«  S)  tkH  y=[2.0t-07,2.0c-07] 
iyNtM(1)+1,0.8*y(1),'BatM  Ma' 
aplat.x,y,calar=3 

;«M  dw  cakAtai  iaaas 

tanat  h  4-IWMA/afe>/1-Pinu,tlwrt  art  aik  8  rtltrme$  kts 

;iivei  la  dw  Ma  far  tfea  targat,  aa  Oaa't  plat  aariaal 
I  (tapgat  aa  S)  tkai  $ 

aplataariaal  ^aariaal  .8aaa(6),payH = 1 
I  (tanat  a«  4)  ap  (tanat  ap  S)  thaa  S 
aiiat,aariaa3j,aariaa3.4aaa(8),payai=;2  $ 
aiaa  aplat,aariaa2.x,aapiaa2.iaaa(8),payai=2 
aplat.aapiea4.z,aariaa4.4aaa(8).payai = 4 
aiiat,aeriaa8.z,aariaa8.daaa(6),iayM= 5 
aai 

2:  laili 

;T1da  aactiaa  ioaa  a  plat  af  aaany  iapaatlaa  varaaa  Z  far  al  rad 
Z  Uvaa  k  caatiawtara 

I 

;alaMr  a  chalca  al  aotpat  prapklca  iavicaa 
pktclMica=iinaaaB(S 

t'Chaaaa  a  praphica  iavica  -'.Xwkdowa',  $ 
'HPGL'.'Pa'.'Paatacrktl.TIdasO) 
caaa  pktcbalca  af 
1:  aat_ptot,Y 
2:  aat  jdat'IT 
8:  aat  ptot,?!!' 

4:  aat  jrtat.'PS' 
aadcaaa 

jlet  tka  Bataa  iata  far  rafaraaca 

I  (tanat  aa  S)  tlwa  bagk 
;tkara  ara  aoly  8  rafaraaca  km 
;plvM  k  iata  far  tkk  tarpat 
pkt  k,raf2.zpipht.ral2.Pasa,d0e=8trtrki(dde),S 
ytik^ltosa  (MaV/kp)M 
idtkslanpRakaal  Dktaaca  (cai)',! 
yraapa=[1.0a-08,1.0a+02],  S 
iraapa=  [0,am(raf2.iripM)],payai=8 
xyaats,raf2.zripM(ralpt8-1 )  1 .8,8 
ral2.Paaa(ralpts-1  ),'r = 2.0cai' 
aaN  aka  kapk 

pkt  k,nlljripH,ral1.«asa,dda=atrtrki(dda},$ 
ytitk='Daaa  (MaV/kp)',t 
xddaslaapkadm  Dktaaca  (cai)',l 
ypaapa=[.00000Q1,aw(nl1.«aaa)],  I 
xraapa = [0,am(raf1  .irlaM)],payai = 8 
xya8ts,raf1  .zriaht(rafpta-1 ) + .08,1 
rail  .dasa(ralpta-1  ),>  s  O.Ocai' 
apkt,raf2.zri8M,ral2.Paaa,payai = 8 
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I  (ekct.neny  m  391.0)  tkM  hegfei 
•lriot,acceirtx1  .acccirtyl  .ktstyli = 2 
•lriot,acceptx2,Mc*pty2,iwstyli = 2 

mM 

I  (tanst  Ml  4)  ar  (tariat  aa  5)  tbaa  I 

xyoiitt,rtl2.iright(rtl|rts-1 ) + .05,1 
ptf2.iasa(retpts-1).'r=2.0cai'  I 

aiaa  $ 

xyaats,rtl2.zriBM(ptr|rtt-1 ) + .05,1 
ral2.iaaa(palptt-1  ),'p = 1  .Ocai' 

aaialM 

apiat.ptl3.2ri|M,pal3.4aaa,payM = 6 
xyMrts,par3.zpiBM(pelpta-1 ) + .05,pal3.iaaa(pa94s-1  ),| 
>=5.0cai' 

aplBt,ral4.zpieM.pal4.iaaa,MyM = 8 
xyMrtt,Ptf4.iPiBM(palRta-1) + .05,ptM.4aaa(ptfptt-1  ),$ 
>s20.0cai' 
x= [18,18] 
y=[1.0a-07,1.0c07] 

I  (target  ag  5)  tkca  y=[2.0c-07,2.0e-07] 
xyaats,x(1)-i-2,0.8*y(1),'Bates  lata' 
aptot,x,y,gsyni=8  ;,color=3 
aglet  [x(1)-10,x(1)],(0.5*y(1),0.5*y(1)] 
xyouts,x(1)+2.0.45*y(1),'PreseBt  iniTRAN  Modaf 
I  (target  ae  5)  aai  (elect_eaergy  eg  391.0)  ttea  hegM 
aglotaccegtx3,accegty3,iMatyie = 2 
aglotacceptx4,accegty4,lDestyk = 2 
eplot[x(l)10,x(l)l,[0.25‘y(1),0.2S*y(1)],| 
beatyii=2 

xyoiits,x(1)+2,0.2*y(1),'Batea  MXQ>T  Modar 

aidH 

jht  the  cafcaiated  dosea 

target  b  4-PMMA/ar/1-PMMA,tkere  are  aaly  3  relereace  bwa 
;glvea  k  tke  data  lar  the  target  ea  daa’t  plet  aerieal 
I  (target  ae  5)  tkea  3 

agletseriea1.irigM,seriea1.dese(6)  ;,payn=1 
I  (target  eg  4)  ar  (target  eg  5)  thei  I 
aglotserieaS.iPigfetaeriea3.deae(8)  I 
eke  egktaerka2.iPigM,aerks2.doae(8) 
egktserka4.irightserkf4.dese(6) 
egkttepks6.zrigM,serka8.deae(e) 


;nk  aectka  dees  a  pkt  el  tke  eoergy  degesUes  cafcakted  at  a 
;dNaea  radhs.  Tlw  tetal  eoergy  degesKed  k  gkned  ever  dw  vabea 
;fer  degesUaa  ly  eaci  type  el  akctrei  ae.  grhury,  seceadary, 
jaod  lOMCkeas) 

cfeokel  swReBadniVlMt  radha  wl  yai  gktT','0.0cM',  I 
'1 .0cai','2.0cai','5.0cai','1 0.0cai','20.0ai',$ 
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'Exr],Tltl(=0) 
csM  dMictl  al 

1:  heik  ftsartasstariaal  ftnMfen=0.0  ftaii 
2:  Itgli  tsarlu=sariBs2  ipaABSs1.0  ftaii 
8:  balk  ftsariatstapiasS  ftpa<n=2.Q  ftaii 
4:  bafki  iaariaa=aariaa1  tnMn=S.O  ftaii 
5:  bubi  ftsariaasaariasS  ftii«M=10.0  ftaii 
6:  baibi  ftaariaa=aariaa8  ftrikKs20.0  ftui 


;alaw  a  cbaica  af  aitpit  imblct  iavicaa 
liatcbaica= wanes 

['Cbaasa  a  rublci  iavica  -'.XiMawa',  $ 
'HPa'.Ta'.'Paatscrkn.TWa = 0) 

iistckifeici  ftf 

1:  aat  piat,X 
2:  aafplat.'HP' 

8:  ufplat.'Pa' 

4:  aarHat.'PS' 

aiicasa 

raiiuB = atrbii  (paibu) 

■axval=MX(sarias.iasa)  ftBbnal=Hbi(tariat.iaaa) 
piat_la,aaries.z,tapias.iasa(0),$ 

tUa = atrtrbn(title) ->■ 'ICRaibis = ' + atPi4i(paiKsi.* ,  J) -i- $ 
•ai'.f 

yUUesDaaa  (MaV/baV.S 
xUUaslaHinaul  iiiat.ii«a  (|/ai2)',S 
yraiia = [abnnlBvrfla.S 
xraiia = [0,6a],]i«aP|k  s  [4.4] 
xyaits,tariap(5)J+2,aariaa(5).iaaa(0),TrbMnr 
aiiat,a?*‘'a«  ^sapias.iasa(2) 
iyaita,aM;'(5).l-»-2.atr'i(5).iasa(2),Kiack-ais' 
aplat,t  ./,aapiaa.iasa(4) 

Kyaats,!:  ai(5).i-«-2,aarias(9).iaaa(4),’G-aac' 

ailat,sartoa.i,8apiaa.iaaa(6) 

xyaits,sariea(5).z+2,sapies(5).iaw(8),latar 


baak 

;ilaw  «  cMka  af  autpat  •rajrtika  iavicaa 
pktcbaica=vnMMi(S 

[’Cbaasa  a  rapbks  icvka  -'.Xwkiawa',  $ 

'HPa.','Pa','Pastscrkt'l,Tlk=0) 

casa  pktcbaica  al 

1:  aat_pkt.T 
2:  aat  pkt.'HP' 

8:  satIbkt.TGl’ 

4:  aat_pfct,Tr 

aiicasa 

;alaiv  isar  ta  Mka  a  cistak  pkl  wU  pracaiiai 

;iata  flk 

stai 
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■mIcsm 

K  (ilitcMct  •«  1)  tbM  i=1 

;!  the  hcvicf  is  a  hariespy  ievica,  ihii  tha  spUaa  si  prihachu 
MMFdcapy 

I  (pistchaics  sa  2)  ar  (pistchaica  aa  3)  ar  S 
(pistchaiea  ea  4)  thea  I 
i=2 
caas  i  af 

1:  pPiBt.'COIII!HUMG  N  X  WMDOWS  MODE' 

2:  ham 

hari;  aatjPat.'x' 

harPeapy = wawaairHaricapy  OptisM',$ 

TpsPica  a  HarPeapy  aa  tha  Ptottar',! 
maha  aa  aPPUsa  ta  tha  mt'.Ha  HapPcapyl.$ 

tXls=0) 

caaa  harPeapy  af 

1:hem 

Mt_piot.'hp‘ 

Pevi^/closa 

lls='iPL‘  +  strlawcaaa(!P.aaBW) 
carf  =  'espy '  -I-  ns  +  '  esaO:' 
apaam,  caiP 
aat.ptot'x' 

aoP 

2:  ham 

priat.  '***  Eater  aay  aPPitisaal  mt  aaaatatieasS 
kterKtively.*' 

prtot  '***  Type  .caatiaae  ta  caathuw  the  prm^.  ***' 
step 

pats,  harp 


sat_mt.'hp' 

Pavica,/cleaa 


prtot,'****  Type  .caathaie  ta  hah  the  prapraai  * 
prtot,'****  ar  aster  istaractiva  WL  caaaaaaPt. 

step 

sstjPetT 

chaica2=Mnaeaa(['Maks  a  chalcar,l 

'Maha  aaather  mt  aPth  carrsat  Pata',  $ 
'Ratara  ta  auto  awas','Exr],T1UB=0) 
cats  chaics2  af 

1:  pata,  hap 
2:  ratara 
aaPcaaa 
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PM 

COMMON  PLOT 

MMM,nMS,nMias,zfact,ai|aii.striM,s«^  ,MriM2,sfritt3,$ 
Mrin4,ttriMS,tarine.stPic2.ilKt_MMriy,tapfit,ptl,ptl1  ,nI2.S 
ptl3,NM.P0Mt>.acMptx1  .accMtyl  ,iicMftx2.acctMy2,$ 
acc»ptt8,»CCTity3,iCCTptt4,icctfty4,z8tirt,|itciiilct,cliilci,| 

tamtsWMMaiNkat  h  tH  TamtT  $ 

•Umtimf,  im  Mimi  by  MmikMi'  .1 

'PMKU  aicli)',S 

'1PMNU,  Mr  fim,  4-PMMA','4-PMMit,  Mr  G«,  1-PMMIll,  TUt^O) 

cast  tapfct  •! 

1:  bdbi  ;far  MnmbiMi  tarftt 

r«l=  {NlM'ncc,Effiarr(5),zriMrt:fliarp(5).raAK:0.0,l 
iMc:niarr(5)  } 

r«l1=rtf  ftrt(2=rtl  Ij<tl3=rtf  ftnM=rtf 
Nlpts=S 

acce|ttx1=ntarr(r«tpt<)  Naccepty1=niarr(rtl|rts) 
Kceirtx2=niarr(rc(pts)  ftMcapty2=fllarr(r«fptt) 
accc(rtx3=ttarr(rtfRts)  ftacc8|rty3=ftarr(r«l|rts) 
aec««tx4=llarr(rcipt*)  NMcepty4=flltpr(rtt|rts) 

Nil  .1= [11 .50,24.50,38.00,51 .50,08.50] 

Nl1.raAis=0.0 

Nil  .zrigtat = [4.58,8.88,1 4.08,20.28,20.14] 

Nl2.IsNf1.Z 
Nl2.ir|ibt  s  Nil  .zrigH 
Nl2.Ndbn=1.0 
Nl3.z=Nflj: 

Nl3.zri9bt = Nil  .zrigM 
Nl3j‘atttt=5.0 
Nl4.Z=Nl1.Z 
Nl4.zrieH = Nil  .zrigbt 
Nf4.raMn=20.0 

I  (diet  Mem  ee  381.0)  Oh  begbi 
Nil  .dOM = [1 8.0,8.5,2.4,1 .2,0.8] 
acccptx1=Nl1.zrlcH 
KCtptyl  =  [20.0,5.4,2.1,1.1,0.85] 
Nl2.den=[3.2e02,1.4e-01,2.3c-01,2.2c-01,1.5e-01] 
accciRx2=Nl2.zri|bt 
acce|rty2= [2.5c-02,.18..82,0.2S.0.14] 

Nl3.iOM =[1.1  e-04,3.2c-04,8.5t-04,2.0e-03,2.5e-03] 
acccptx3=Nl3jri|bt 

acceety3= [4.0c-04,2.0c-04,8.0e-04,1.8e-03,1.8e-831 
Nl4.deM= [1 .8e-08.3.2e-06,5.0e-08,B.0e-08^8.0e-0e] 
acceptx4=Nl4.zriMM 

accepty4 = [5.4e-0e,1 .8e-07,8.5e-08,5.5e-08,5.0e-0e] 


I  (iiect_Nam  >4  200.0)  tbN  beibi 
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Ptd  .riMt= [S.0,1 .8.075,0.35.0.181 
Kl2.iM«= [0.58,0.24,0.18,0.1.0.055] 
l>•l3.iM•= [1 .4t-04,4.3a-04,1  .lt-03.1 .8t-03,2.05t-03] 
nM.8mc= [3.31-06,1 .1  ■•05,1 .3^05,1 .85i-05] 


2:  belli  ;tar  Lui  fiimitil  by  Umkmt 

pels  {pcfM‘«Bcc,Elllapp(5),zpi[8tt:flUfT(5).nAn:0.8,$ 
i■s■:llapp(5)  } 

p«n=Ptl  apcl2=pcl  lPtl3=Nl  ir«f4=PCf 
kI|U=5 

icci|lx1=atapp(ptlpts)  &icc«|ty1=IKiPP(Nl|U) 
KM|tx2=HUpp(p«llt<)  taccapty2=fltipp(pei^t) 
acmMx3=fltapp(r«fpts)  iacecity3=llipp(pa0ttt) 
■cci|tx4=l8ipp(i>al|ts)  ldKCipty4slltapp(pflMt) 

Kf1  [10.8,25.0,38.0,51 .0,67.1] 

Nil  .ZPilbt = [1 .28,8.81,1 1 .83,17.05,23.32] 

Nn.PaAls=0.0 

Nl2.z=Ptn.z 

pel2.ziH|bt = Ptfl  .zpiilit 

P«l2.PldbUs:1.0 

p«(3.z=ptl1.z 

Nl3.ZPiBht = Pill  .ZPilbt 

Ptt3j>adbis=5.0 

P«I4.Z=P«I1.Z 

pcl4.zpi|bt=pcf1  .ZPilbt 

pcM.Pidbis=20.0 

I  (■ket  eiepnr  d  381.0)  tim  beibt 
P«I1  .im=  [21 .0,3.2,1 .3,0.48,074] 
icctptxl  =p«(1. ZPilbt 
iCMityl  =  [30.0,2.85,2.0,0.55.075] 
pcl2.l«ns  [0.18,0.21,0.10.0.14,8.0t-02] 
icci|tx2=Pfl2.zpi|bt 
icc«pt](2 = [0.2,0.3,074,0.1 3,0.1 ] 

Ptl3.l«M= [4.5t-04,1.7»-03.2.5i-03,2.8e-03,2.8t-03] 
ac€C|tx3=Ptl3.ZPi|bt 

tcc«|ty3= [3.8e-04,2.2«-03.4.1»-03,5.5e-03,3.0e-03] 
Pil4.4Me= [5.0c-06.1.5«-05.2.0e-05.27*-05,2.1t-05] 
icc«|tx4=p«f4jpi|bt 

aCM|ty4= [4.2»-06.1 .3e-05.1 .5t-05,7.8«-08,2.8c-08] 

I  (■het  MWPiy  §4  200.0)  dwi  b«|bi 
Ptt1.ioM=  [7.5,.80..54..34,.41] 

Pil2.iaM= [.28,71, .10,.055,.03] 

Ptl3.lm= [5.8C-04.1 .8t-03,2.0i-03.2.4f-03,27i-03] 
Pff4.isu= [8.6e-08,1.7i-0S,2.3c-05,2.5t-05,2.8l-05] 


3:  bcibi  ,1ap  PHnu  Mriy 
pels  {ptt»PMCi,z:H«pp(5),zpbkt:flt«pp(5),pidbw:0.0,| 
ieti:fltipp(5)  } 
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rfllspti  trtlZ=rtl  trtl3=rtl  lirtf4=rtf 

ptIVttsS 

acctinx1s;Rarr(ptlint)  ftacc«ity1=Sarr(rtl|dt) 

acca|rtx2=niarr(rtlptt)  ftKciityZ=ltarr(ptf|ttt) 

iccairtx3=niarr(r«lpts)  iacctpty3=:Rarr(p«lptt) 

acca^4=Rarr(ptlpts)  iacc«pty4=Rarr(rtlpts) 

rtn  .z= 110.0.20.0.30.5,41 .0.51 .0] 

rtn  .zri|M= [8.72.17.44.20.57.35.71,45.10] 

rtn.rMlw=0.0 

Ptl2.z=rtl1.z 

rtl2.zr|Rt=rtn.zri|ht 

ptl2.ridin=1.0 

ptl3j=rtnj 

rtl3.zri|M = rtn  .zriitt 

N(3.ri4i>s=5.0 

ptl4.z=rtl1.z 

ptM.zri|H=ptfl.zrieU 

rtM.railBt=20.0 

I  (elect  eeeny  ee  301.0)  tkei  heih 
refl  .im = [1 0.0.4.8.1 .7,0.75,0.21  ] 
acceitx1=ref1.zri|M 
eeceptyl  -  [1 9.0,2.0.1 .5,0.8..38] 
Ne.«eee=[3.2e-02,1.4e-01,2.3e-01,2.2e-01,1.5e-011 
ecceptx2 = re(2.zri|M 

eccepty2= [5.8e-03.8.0e-02,0.28,0.30,0.70] 

rel3.4ese=[5.0e-05,1.8e-04,4.8e-04,1.1e-03,1.7e-03] 

ecce|rtx3=rel3.zriRrt 

eccepty3s  [2.0e-05,5.4e-05,2.0e-83.2.4e-03.6.0e-03] 
peM.deee = (2.0e-08,3.8e*06,5.1  e-06,8.5e-08,7.1e-081 
ecceptx4=pef4.zri|M 

ecceeiy4 = [2.4e-08,1  .Oe-08,1 .0e-08.e.0e-07.3.1e-06] 


I  (elect  eeeriy  ee  200.0)  dm  talk 
red  .itu= [8.0,2.0,.85,.92..88] 
pe(2.tasc= [.059,.22,.16..09..05] 

Nf3.taee= [1 .1  e-04.4.4e-04,1 .3e-03,3.2e-03,4.9e-031 
reR.taee= [1 .2e-05.1 .5C-05.1 .7e-05,1 .8e-05,2.0e-051 


4:  talk  ;1PMIU  lelewel  ly  m  *  W4  kleweO  ky 

M-PMMA 

rel=  {rclercicc,z:nterr(8),zrl|M:IRiiT(8),ptdke:0.0,8 
leec:flterr(8)  } 

ptflsrel  lrel2=pc(  trel3=pel  &Kl4=pe( 
refitesB 

iccc|tx1=IKirr  (rente)  8kcce|ty1=Rirr  (rente) 
acce|rtx2=flterr(renti)  fteccepty2=IKirr  (rente) 
ecceitx3=Rirr  (rente)  ieccepty3=Hterr  (rente) 
icce|rtx4=fltirr (rente)  tecceptylsllirrtrente) 
prkt.nrmts'CEiter  the  kiitb  k  the  elr  m  k  I 
ca  -M)' 
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Nn.l=  [10.35.10.40,20.90,31  .ea42.3a53.101 
Ptn.iri|ht= [0.72,48.02,50.74, e5.4a74.ia82.90] 
rtllJ‘t0ks=0.0 

PtII  .«nt= [18.a0.55,0.5a0.45,0.3a0.25] 

acciptxisrtlijrilfet 

accwtyl = [18.a,2a.4a-32..S7,.281 

rtl2.isNl1.z 

rfl2.zri|M  3  Nil  Jrfitt 

Nl2j‘iMns2.0 

Nl2.iM*= [1.8«-0ai.7»-02,2.8t-0a4.5c-02,9.0>42,4.5c-02] 
acnptx23Nl2.iri|ht 

icc«9ty2= [8.2l-04,2.8t-82,8.0t-02,0.5t-02,7.0t-02,a5l-02] 
Nl3.ZsNn.Z 
Nl3.zriiM = Nn  jriitt 
Nl3.NMn=5.0 

Nl3.dNl3  [2.8c-0a8.0l-04,8.5t-04,1 .3i-0a2.5c-0a4.0t-031 
ac€tptx3=Nl3.zri|H 

acctptyS = [4.0c-05,2.0c-0a2.2«-03,3.8«-0a7.5l-0a0.0i-03] 
nI4.Z3NI1j 
Nl4.ZPi|M = Nil  .Zrigtt 
Nl4.NdiM=20.0 

Nf4.iMt= [4.5«-07,4.0c-05,7.0c-08,8.0c-0ai  .0l-0ai-5t-051 
aeetptx4=Nl4.zriiH 

KCipty43  [2.1  •-07,5.5t-05.1 .04-00,1 .3t-0a2.0c-0aa54-0e] 

4Ml 

5:  bilk  ^4-PMMA  fiOoinO  by  M  kr  lai^  ItliwiO  br  1-PIMU 
Nf  3  {NliN4C4.z:IOarf(8).zri|bt:fli4rr(fl),NkwJ.at 
itM:marr(e)  } 

NfIsNf  iNl23Nf  tNl33Nf  tNl43Nl 
NlbtS38 

prbai4nnt3'(‘EMtp  Um  Imtb  •!  tbt  *  up  k  0 

m  -M)' 


Nil  .z  3  [1  o.a2o.a3i  .a4i  .a43.i,53.i  i 
Nn.ZPl|bt3  [8.72,17.44,28.ia34.88  74.ia82.90] 
NflJ>a0kS3  0.0 

Nil  J4M3  [o.oao.oao.oao.oao.oao.oo] 

Nl2j3Nt1J! 
nI2  jpWm  3  Nil  Jrikn 
Nl2.N4bR3  2.0 

nI2  J4lt3  [1.8t-03,5.8t-0ai  .84  02.8.0t-02,2.l4-0a3.l4-02] 

Nf8j3Nf1j 

Nf8Jri|bt3Nfljrkbt 

Nl3j>a*S35.0 

nI8  Jm3  [8.0t-0aa74-04.e.4t-04.1 .54-0ai.04-02.1.1t-02] 
nI4J3nI1j 
Nf4  jrkbt  s  Nil  JPkbt 
Nl4j'a«n320.0 

Nl4.i4S43  [2.84-0a7.04-0e.8.1t-Oa8.9t-Oaa9t-O4,O.94-05] 
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APPENDIX  C.  SUPPLEMENTARY  FIGURES 


The  supplementary  figures  are  intended  to  enhance  the  discussion  in 
chapter  3,  but  are  not  essential  to  understanding  the  results  of  this 
work.  The  figures  begin  on  the  following  page. 
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( \OB\a  1 6>i  /  AeN  ]  esoQ 


Figure  Cl.  Dose  Breakdown  on  Beam  Axis  for  Target  of  Figure  11  with 
10000  Electrons.  On  the  centerline  the  dose  is  almost  completely  due  to 
primary  electrons.  Note  the  lack  of  dose  after  a  radiation  length,  which 
is  34.4  cm  in  PMMA. 


(  pe|9  /  &>|  /  A9IAI  }  esoQ 


Figure  C2 .  Dose  Breakdown  at  2cin  Radius  for  Target  of  Figure  11  with 
10,000  Electrons.  There  is  very  little  energy  deposited  in  the  first 
20cin.  Note  the  reduction  in  scale  from  the  previous  figure,  i.e  the  dose 
is  reduced  guickly  off-axis. 
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Dose  Breakdown  -  r=5cm,  4-PMMA,  39cm  Air,  1 -PMMA 


( IDeie  /  6^  /  A9lAj )  esoQ 


Figure  C3.  Dose  Breakdown  at  Scm  Radius  for  Target  of  Figure  11  with 
10,000  Electrons.  Dose  is  negligible  until  electrons  reach  the  air  gap. 
Photon  processes  obviously  play  a  more  important  role  off-axis. 


Longitudinal  Distonce  (cm) 


OoM  from  Primqrv  Ej^otrom  ~  fZOcm.  4>  PliM.  3^cm  Air,  t  ■  PIbftMk 


Figure  C4.  Dose  Breakdown  at  20cTn  Radius  for  Target  of  Figure  11  with 
10,000  Electrons.  Note  reduction  in  scale  from  previous  figures.  Also 
note  that  the  dose  from  primary,  gamma-secondary,  and  knock-on  electrons 
is  roughly  equal. 
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(  P9|9  /  6>|  /  ASIA]  )  aSOQ 


Figure  C5.  Low  Resolution  391  MeV  Dose  Distribution  for  Target  of  Figure 
4.  Comparison  of  CYLTRAN  calculations  using  2  different  random  number 
seeds,  all  other  parameters  remaining  the  same. 
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Figure  C6.  High  Resolution  Dose  Distribution  for  Target  of  Figure  10 
with  Different  Random  Number  Seed.  Compare  to  figure  10  to  see  that  the 
random  number  seed  can  cause  the  results  to  fluctuate  about  20-50%. 
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RwoKition  Mod<l 


Soan  Of  A(r  Go 


OoM  from  Kno<t<-«n  OaciiWM  -  r-20cm,  Wah  WteJulion  M<mW 


Figure  C7.  Dose  Breakdown  at  20cin  Radius  for  Target  of  Figure  10  with 
Different  Random  Number  Seed.  Notice  the  small  number  of  collisions 
experienced  in  the  air  gap.  The  change  in  random  number  seed  causes  a 
drastic  change  in  the  distribution  in  the  air. 
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